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Preface 


On June 28-30, 1993, the Workshop on Atmospheric Transport on Mars 
was held in Corvallis, Oregon. The workshop was organized under the 
auspices of the MS ATT (Mars Surface and Atmosphere Through Time) 
Program of NASA, and was jointly sponsored by the Lunar and Planetary 
Institute, Oregon State University, and the Oregon Space Grant Consortium. 
More than 50 scientists (including a number from Europe) attended the 
workshop, which was the first such meeting to focus upon circulation pro- 
cesses in the Mars atmosphere. The timing of the workshop placed it almost 
on the eve of the arrival of Mars Observer at Mars, so that the presented 
papers gave a picture of the “state of the art” in Mars atmospheric science 
just prior to the expected arrival of new data. (Unfortunately, the loss of Mars 
Observer has now pushed the expected acquisition of such data at least 
several years into the future.) The workshop highlighted a host of recent 
advances in atmospheric modeling and analysis — advances that will be rel- 
evant to any future observations. 
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Summary of Technical Sessions 


SESSION I: ATMOSPHERIC DYNAMICS 
AND CIRCULATION 

It is, of course, winds that act to produce atmospheric 
transports of various constituents such as dust and water and 
trace-chemical species. Carbon dioxidemakes up most (95%) 
of the mass of the atmosphere so that, essentially, its transport 
is the wind (density weighted). As the winds redistribute the 
total atmospheric mass on the planet, the surface-pressure 
field reflects these changes. The winds are fundamentally 
driven by the differential solar heating of the planet, which in 
turn can be affected dramatically by the distribution of atmo- 
spheric dust There is also, of course, a contribution to the 
winds associated with the condensation and sublimation of a 
sizable fraction of the atmosphere in the polar regions. Un- 
derstanding the winds (that is, the atmospheric circulation 
and its underlying dynamics) that are responsible for trans- 
ports was the focus of the opening session of the workshop. 

The session opened with an invited paper by C. Leovy on 
Mars atmospheric circulation. Leovy focused his discussion 
on three aspects of the circulation: heat transport into the 
polar regions, thermal tides, and gravity waves. He began by 
emphasizing the extent to which recentmodeling efforts have 
outstripped the available observations, highlighting the criti- 
cal importance of Mars Observer-like data, hi relation to the 
heat transports, Leovy focused upon the north-south asym- 
metry that general circulation modeling experiments have 
revealed. This is associated with the fact that the atmosphere 
is dustier near the time of northern winter solstice (as well as 
the fact that this occurs near perihelion), and with the model 
result that the transient eddies are considerably weaker dur- 
ing southern winter than northern winter. The latter appears 
to be at least partly a consequence of topographic differences 
between the hemispheres. The topography is poorly known at 
present, however, so the model results must be regarded 
skeptically; obtaining good topographic data is vital for in- 
vestigations of the atmospheric circulation. In relation to 
thermal tides and gravity waves, Leovy emphasized that both 
may be dominant players in the middle atmospheric region. 
From a modeling perspective, gravity waves present a diffi- 
cult challenge since they cannot be resolved without very high 
resolutions and the current parameterizations of their effects 
are rather crude. 

Following Leovy ’s review talk, there were four papers that 
dealt with observations of the atmospheric structure and 
circulation. M. Santee discussed analyses of the thermal 
structure and circulation based upon Mariner 9 IRIS data for 
late northern winter. One of the most interesting aspects of 
these results was that they indicate that the north polar region 
was warmer at high levels (above about 40 km) than both low 
latitudes and the south polar region. This implies a diabatic 


mean circulation characterized by sinking motions at high 
northern latitudes and rising motions in low latitudes. The 
existence of such a mean meridianal circulation implies 
strong forcing by eddy motions — on the order of 50 m/s/day 
in the winter polar region. Unfortunately, the Mariner 9 data 
do not indicate what particular type of eddy motions may be 
responsible for this forcing. W. Mc Millan then discussed the 
retrieval of atmospheric temperature profiles from Mariner 9 
IRIS data as part of an effort to develop and test a fast and 
efficient inversion scheme for use with Mars Observer TES 
spectra. L. Martin gave a brief review of observations of 
clouds and cloud motions, obtained telescopically and from 
spacecraft. To the extent that such motions can be imaged by 
spacecraft, they can be very valuable as direct wind indica- 
tors, but the interpretation of such motions is often ambigu- 
ous. Finally, E. Chassefi&re gave a talk on some of the results 
from the Phobos 2 solar occultation measurements. These 
results, for cloud, dust, and ozone abundances in the middle 
atmosphere, provide constraints on the strength of the verti- 
cal mixing in this region. 

Following the observational papers, the remainder of the 
papers discussed recent results from a large variety of mod- 
eling studies of the martian atmosphere. Whereas as recently 
as two years ago the NASA/Ames Mars General Circulation 
Model (GCM) was the only such model for another planetary 
atmosphere, several such models are now under develop- 
ment In particular, there is now a NASA/GISS GCM being 
developed (reported on by M. Allison), as well as a NOAA / 
GFDL GCM (reported on by R. J. Wilson). Additionally, 
there is a French (LMD) GCM effort and a British GCM 
project (these were reported on by F. Hourdin and S . Lewis in 
later sessions). There are also several “simplified GCM” 
studies underway, as reportedon by R. Haberle andM. Collins. 
The study discussed by Haberle was particularly interesting 
in that continuous multiannual (10 Mars years) simulations 
had been performed. These appear to show sizable interannual 
variations — in the absence of topography and interactive dust 
transport — that peak in the southern summer season . [Haberle 
also presented some intriguing (but baffling) results on the 
apparent correlations between global dust storms on Mars 
and perfect games pitched in major league baseball.] The 
study discussed by Collins focused upon the transient eddies, 
and the (sometimes) “regular” nature of their structure. The 
NASA/GISS GCM is interesting in that it includes param- 
eterized gravity wave effects, and Allison discussed the na- 
ture of the zonal wind balances in the model. The NOAA/ 
GFDL GCM has been used primarily for experiments focus- 
ing upon thermal tides and the nearly resonant Kelvin wave 
modes. 

The remainder of the modeling papers discussed more 
specialized and/or nonglobal circulation models. Chassefifcre 
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presented results from simulations of the middle atmosphere 
circulation with a zonally symmetric model incorporating 
gravity wave drag. J. Hollingsworth then discussed linear 
modeling of forced stationary waves in the winter hemi- 
spheres of Mars. M. Joshi reported on the examination of such 
waves, having a boundary current character, in simulations 
performed with the British GCM. A model of the circulation 
in the very high reaches of the atmosphere, the thermosphere, 
was discussed by S. Bougher. It is hoped that this model can 
be noninteractively coupled with the NASA/Ames GCM to 
yield simulations of the atmosphere from the ground well up 
into the thermospheric region. Finally, T. Siili and T. Parish 
reported on results from mesoscale modeling studies. The 
former study has examined circulation driven by gradients in 
surface properties, while the latter has focused on katabatic 
flows in the north polar region. 


SESSION E: DUST 

Dust can have a tremendous influence upon the heating of 
the martian atmosphere, which in turn drives the winds that 
transport and raise the dust Understanding the raising, trans- 
port, and deposition of dust is thus crucial to a full under- 
standing of the atmospheric circulation and climate. J. Murphy 
began the session on dust with an invited paper discussing 
dust transport He began by outlining some of the key ques- 
tions relating to the dust cycle, and briefly reviewing some of 
the available observations that bear upon dust transports. 
Murphy then discussed some of the results from simulations 
of dust transport and lifting that he has been carrying out with 
the NAS A/ Ames GCM coupled to an aerosol transportmodel. 
The simulations in which lifting is allowed to take place are 
especially interesting, despite the uncertainties in the quan- 
titative nature of the lifting process, and appear to show both 
positive and negative feedback. The simulations evidence 
considerable lifting in the winter hemisphere in association 
with midlatitude transient eddies, and appear to show that 
topography acts to reduce the total amount of dust raising 
(because of the presence of sizable regions at relatively low 
pressures). 

Following Murphy’ s review were four papers dealing with 
various observations of dust First, P. James discussed Viking 
imaging observations of dust activity in the northern hemi- 
sphere during spring. There was one observation of a quite 
extensive (regional-scale) storm at high latitudes in the spring 
season following the 1977 fall and winter global dust storms; 
the frequency and size of such storms could be important for 
the global dust cycle. T. Martin presented dust opacity maps 
derived from Viking IRTM data. S . Lee and R. T. Clancy then 
reported on dust particle sizes and compositions and studies 
of dust deposition and removal from the surface. Improving 
our knowledge of the latter is essential to gaining a better 
understanding of dust raising and transport. 


SESSION IE: VOLATILES 

The C0 2 and water cycles are intertwined with the dust 
cycle, and these were the focus of the next session in the 
workshop. B. Jakosky led off the session with an invited 
review of these two cycles. He discussed the fact that both the 
dust and water cycles evidence very large interannual vari- 
ability, but for the most part the C0 2 cycle does not The 
observed changes in the shape and structure of the retreating 
seasonal polar caps are evidence of interannual variations in 
the C0 2 cycle. Though the existence of one permanent C0 2 
cap seems understandable (as opposed to two), the funda- 
mental reasons why that cap exists in the south at present are 
not at all clear. Finally, the longer-term behavior of the polar 
caps is not well understood at all. In the case of the water 
cycle, Jakosky emphasized the very large interannual vari- 
ability that is evidenced by Earth-based observations, and the 
basic question of whether or not a net transport from north to 
south is currently occurring. 

Following the review, there were two papers dealing with 
recent GCM simulations of the annual pressure variations, 
given by J. Pollack (the NASA/Ames GCM) and Hourdin 
(the LMD GCM). Both studies have shown that there is a 
sizable “dynamic” component of the seasonal variation in 
surface pressure on Mars. This component is associated with 
the redistribution of atmospheric mass, in conjunction with 
seasonal changes in the zonal-mean circulation. Pollack 
showed that the atmospheric heat transport is very important 
in producing the observed strong asymmetry in the annual 
pressure curves, and this is largely due to greater atmospheric 
dust loading during northern winter. This enhanced dusti- 
ness may play arole in the formation of C0 2 ice clouds in the 
north polar region in winter, evidence (from Viking IRTM 
temperature data) for which was then discussed by F. Forget. 
A. Zent followed with apresentation of somerecentmodeling 
results for the transport of water in the atmospheric boundary 
layer. Finally, Chassefifcre and J. Barnes both touched upon 
the topic of eddy mixing. Chassefifcre reported on a photo- 
chemical modeling study that suggests a lower limit for the 
strength of the vertical mixing, and Barnes discussedaproject 
to characterize the nature of tracer transport — both diffusive 
and advective — that is implied by the Ames GCM-simulated 
atmospheric circulation. 

SESSION IV: MARS OBSERVER AND FUTURE 
SPACECRAFT MISSIONS 

The papers presented during the first two days of the 
workshop emphasized the fact that the wide range of models 
that have been developed have outstripped the available data. 
The final session of the workshop dealt with Mars Observer, 
as well as with other possible future spacecraft missions. 
Despite the loss of Mars Observer, the analysis approaches 
discussed are relevant to future global, remotely sensed data- 
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sets. R. Zurek began things with an invited review of the 
various atmospheric measurements planned for Mars Ob- 
server. The PMIRR measurements of temperature, dust, and 
water were perhaps primary, but extremely important obser- 
vations were also to be made by TES, radio science, MOLA, 
and MOC. The opportunities for synergistic studies would 
have been considerable. Four papers that dealt with various 
aspects of planned analyses of Mars Observer data to infer 
winds and transports were then presented. Lewis reported on 
the progress of an effort to develop a GCM to be used for data 
assimilation with PMIRR temperatures. This project uses the 
data assimilation code employed by the British Meteorologi- 
cal Office. D. Banfield then presented results from a study of 
optimal gain functions for data assimilation using a simpli- 
fied dynamic model. A planned effort to produce synoptic 
maps from Mars Observer data using a simplified (linear- 
balance) atmosphericmodel was thendiscussedby H. Houben. 
Finally, a diagnostic scheme designed to infer winds and the 
small-scale (e.g., gravity wave) zonal drag from the PMIRR 
temperature data was discussed by Barnes. 

The workshop concluded with several papers dealing with 
possible post-Mars Observer spacecraft missions that could 
yield advances in our understanding of atmospheric trans- 
port Haberle gave a brief synopsis of recent developments in 
the MESUR mission planning process (in particular, sharp 
reductions in the size of the surface network), and then 
opened the floor to a discussion of the impacts of these 
developments on atmospheric science as well as ways to try 
to retain as much of this science as possible in the MESUR 


mission. P. Smith followed with a report on the IMP instru- 
ment that was recently selected for the MESUR Pathfinder 
mission that is planned to precede MESUR. This imaging 
system has the capability of making a number of measure- 
ments of atmospheric properties. S. Limaye gave a presenta- 
tion on two possible orbiter mission concepts, MOES and 
MOMS. The first could obtain complete local time-of-day 
coverage of Mars, while the second would be a means of 
obtaining temperature measurements without any need for 
sounding instruments (through radio occultations between 
several satellites). B. Clark then presented the final paper of 
the workshop, discussing a possible long-lived balloon mis- 
sion to Mars. 

At the conclusion of the workshop, participants generally 
felt that it had been successful. Because this was the first 
meeting to focus on circulation processes, it brought together 
more specialists (including several from outside the U.S.) in 
Mars atmospheric sciences than had any previous meeting. 
The workshop showcased the recent developments in model- 
ing of Mars global atmospheric circulation. The number of 
global models that now exist is remarkable, but it was clear 
at the workshop that these wonderful tools have gone beyond 
the available data. It was certainly hoped that the next such 
gathering of Mars atmospheric scientists will be one at which 
the data are ahead of the models. While the loss of Mars 
Observer may have dealt those hopes a strong blow, the 
meeting certainly did serve to emphasize the vital importance 
of acquiring global observations of the atmosphere. 
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Abstracts 


THE NASA/GISS MARS GENERAL CIRCULATION 
MODEL: PRELIMINARY EXPERIMENTS. M. Allison*, 
M. A. Chandler 2 , A. D. Del Genio*, A. Lacis*, D. Rind*, W. B. 
Rossow*, L. D. Travis*, and W. Zhou 3 , *NASAGoddardInstitute for 
Space Studies, 2880 Broadway, New York NY 10025, USA, 
University Corporation for Atmospheric Research, NO AA Global 
Change and Climate Program, Institute for Space Studies, New 
York NY 10025, USA, 3 Hughes STX Corporation, Institute for 
Space Studies, New York NY 10025, USA. 

The NASA/GISS Mars General Circulation Model (GCM) is an 
adapted version of the GISS Global Climate/Middle Atmosphere 
Model [1], specifically developed for the diagnostic validation and 
objective analysis of measured atmospheric temperatures from the 
Mars Observer Pressure Modulator Infrared Radiometer (PMIRR) 
experiment The GISS Mars GCM has 23 vertical layers extending 
from the surface to approximately 80 km altitude, representing a 
vertical resolution of about 0.3 scale heights. The primitive (verti- 
cally hydrostatic) equations are solved in finite difference form on 
the Arakawa B grid, with a horizontal resolution of 8° x 10° 
(latitude-longitude). The model includes a diurnal solar cycle, heat 
transport within a two-layer ground, and a high-order “slopes- 
scheme” for the advection of heat in the upper atmosphere [2]. The 
radiative transfer scheme is based on the correlated k distribution 
method for the treatment of nongray gaseous absorption, thermal 
emission, and multiple scattering [3], including options for sus- 
pended dust. A special feature of the model of particular importance 
for Mars is a parameterization of gravity-wave-induced drag incor- 
porating orographic forcing, wind shear, convection, and radiative 
damping. The implementation of the GISS Mars model includes 
global maps of topography, roughness, and albedo. 

In a first series of experiments, excluding phase changes and a 
hydrological cycle, with a pure COj atmosphere, the martian circu- 
lation was simulated for the fixed solstice and equinox seasons. 
The fixed equinox (Ls =0°) simulation includes westerly 40-50 m s~* 
jets at midlatitudes and easterly upper-level flow at the equator 
exceeding 30 m s~*. A hemispheric asymmetry in the zonal wind, 
temperature, and meridional stream function is apparent, with the 
rising branch of the Hadley circulation centered at about 15° to the 
south of the equator and a peak stream function of about 2.5 x 
10 9 kg s-*. The asymmetry is evidently the result of an enhanced 
surface heating in the southern hemisphere associated with its 
relatively elevated topography. For the fixed northern summer 
(Ls = 90°) simulation, the zonal velocities, associated temperature, 
and Hadley circulation (rising in the north and descending in the 
south) are all somewhat stronger, with largely easterly flow in the 
northern hemisphere. For the northern winter (Ls = 270°) simula- 
tion, the zonal velocity and cross-equatorial Hadley circulation 
pattern is reversed, with westerly 70 m s-* flow to the north, peaked 
in the vicinity of the 0.1-mbar level, easterly 90 m s* 1 upper level 
flow to the south, and a mass stream function exceeding 1 x 10*° kg 
s -*. An isolated 20 m s-* westerly jet also appears near the surface 
at this season at southern midlatitudes. A preliminary series of 
simulations including the radiative effects of fixed airborne dust 
were also conducted and generally show still stronger zonal winds, 
sensitively dependent upon the amount and distribution of opacity. 


Following the fixed season experiments, an extended running inte- 
gration of the Mars GCM was initiated, accounting for the variable 
solar declination and distance, beginning approximately 55 sols 
before the southern summer and extending for nearly 400 sols of 
continuous simulated time. The variable season experiments ex- 
hibit a temporal lag of several sols with respect to the fixed season 
simulations, attributed to the thermal inertia of the deep ground 
layer, but otherwise show a qualitatively similar wind and tempera- 
ture structure. 

References: [1] Rind et al. (1988) J. Atmos . Sci, 45, 329. 
[2] Russell and Lcmer (1981) /. Appl Meteor,, 20, 1483. [3] Lacis 
and Oinas (1991) JGR, 96, 9027. 


OPTIMAL GAINS FOR A SINGLE POLAR ORBITING 
SATELLITE. D. Banfield*, A. P. Ingersoll*, C. L. Keppenne 2 , 
* California Institute of Technology, Pasadena C A 91 125, USA, 2 Jet 
Propulsion Laboratory, Pasadena CA 91109, USA. 

Gains are the spatial weighting of an observation in its neighbor- 
hood vs. the local values of a model prediction. They arc the key to 
data assimilation, as they are the direct measure of how the data are 
used to guide the model. As derived in the broad context of data 
assimilation by Kalman [1] and in the context of meteorology, for 
example, by Rutherford [2], the optimal gains are functions of the 
prediction error covariances between the observation and analysis 
points. Kalman [1] introduced a very powerful technique that allows 
one to calculate these optimal gains at the time of each observation. 
Unfortunately, this technique is both computationally expensive 
and often numerically unstable for dynamical systems of the mag- 
nitude of meteorological models [3] , and thus is unsuited for use in 
PMIRR data assimilation. However, the optimal gains as calculated 
by a Kalman filter do reach a steady state for regular observing 
patterns like that of a satellite. In this steady state, the gains are 
constants in time, and thus could conceivably be computed off-line. 
These steady-state Kalman gains (i.c., Wiener gains) would yield 
optimal performance without the computational burden of true 
Kalman filtering. We propose to use this type of constant-in-time 
Wiener gain for the assimilation of data from PMIRR and Mars 
Observer. 

We will discuss the ways in which our technique differs from the 
method of data assimilation commonly used in Earth applications, 
Optimal Interpolation (OI). These differences include assumptions 
usually made in the implementation of OI about the shape of the 
gains, and the way multivariate covariances are handled. Implemen- 
tations of OI typically make the assumption that the gains are 
isotropic in space [c.g., 4]. Propagation of atmospheric waves and 
the westward-marching nature of the observing pattern of Mars 
Observer cause significant asymmetries (about the observation 
point) in the shape of the optimal gains. We therefore retain the 
anisotropic nature of the gains with marked improvement in perfor- 
mance in our models over the isotropic gains of OI. Since PMIRR 
will only observe one dynamical variable (temperature), analyses 
must include cross-correlations between temperature and the other 
variables. A balance (e.g., geostrophic) is often assumed in the 
cross-correlations of many OI implementations [e.g., 4]; however, 
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we impose no explicit balance assumptions on our gains. The actual 
statistical covariances constitute the cross correlations that link the 
dynamical variables in our technique. Thus, the appropriate balance 
between the dynamical variables of the system will be maintained. 

The optimal gains arc, in general, functions of the latitude and 
longitude of the observation point Topographic relief is greater on 
Mars than on Earth, and thus one would expect local topographic 
Gongitudinal) effects to be significant in determining the shape of 
the optimal gain functions. However, the Wiener gains we will use 
in assimilation of PMIRR data will be functions of only the space- 
craft latitude, ignoring their longitudinal dependence. In our 
models, performance is only marginally affected by making this 
simplification as compared to the isotropic assumption discussed 
above. The benefit is that it greatly reduces the number of unique 
gain functions that must be maintained. 

We have developed a method of determining these optimal 
Wiener gains through an iterative procedure. As there are very few 
direct observations of martian weather and none yet that will suffice 
in providing a dataset for use in assimilation of global phenomena, 
we use a randomly forced model as a surrogate for Mars. We then 
assimilate observations of this 'Truth” model into a predictive 
model using a zeroth iterate guess of the gains that are derived from 
the weather covariances of the modeL As the assimilating continues, 
we compute the prediction error covariances that result, as we know 
both the prediction state as well as the 'Truth” state for this set of 
systems. These prediction error covariances will then yield a next 
iterate, which will be a better estimate of the true steady-state 
optimal (Wiener) gains. This technique is then repeated until the 
gains converge to the actualoptimal Wiener gains. In all ourmodels, 
we found convergence in about three iterations of this procedure. 
Furthermore, in comparing with Kalman filters when we could 
implement them, the Wiener gains derived with this technique 
converged to the same functions as the steady-state Kalman gains. 

To back up these statements, we have used a set of simple (and 
not so simple) models tuned to behave like the martian atmosphere. 
The anisotropies in the steady-state Kalman gains (Wiener gains) 
are well demonstrated through a one-dimensional one-variable 
Rossby wave model that simulates a latitude circle on Mars. The 
convergence of the above iterative procedure is also easily seen with 
this model. A one-dimensional three- variable Rossby gravity wave 
model was used to test the multivariate nature of the gains when only 
one dynamical variable is observed. This same model was also used 
to test the effects of topography on the gains. Mars-like topography 
was introduced into the model, and the predictive performance was 
tested with different assumptions about the longitudinal dependen- 
cies of the gains. Finally, a spectrally represented shallow- water 
model on a sphere [5] was used to test this technique with a non- 
linear model of similar (but less) complexity as a GCM. Experi- 
ments with these models suggest to us that this technique should be 
portable to a GCM with no major conceptual changes. 

References: [1] Kalman R. E. (1960) J. Basic Engineering 

(ASME),82D, 35-45. [2] Rutherfordl. D. (1972) JAS, 29, 809-815. 
[3] Bierman G. J. (1977) Factorization Methods for Discrete Se- 
quential Estimation, 244, Academic, New York. [4] Bergman K. H. 
(1979) Monthly Weather Review, 107, 1423-1444. [5] Keppenne 
C. L. (1992) Icarus, 100, 598-607. 


A DIAGNOSTIC MODEL TO ESTIMATE WINDS AND 
SMALL-SCALE DRAG FROM MARS OBSERVER PMIRR 
DATA. J. R. Barnes, College of Oceanic and Atmospheric 
Sciences, Oregon State University, CorvallisOR97331-2209, USA. 

Theoretical and modeling studies indicate that small-scale 
drag due to breaking gravity waves is likely to be of considerable 
importance for the circulation in the middle atmospheric region 
(-40- 100 km altitude) on Mars [1,2]. Recent Earth-based spectro- 
scopic observations have provided evidence for the existence of 
circulation features — in particular, a warm winter polar 
region — associated with gravity wave drag [3]. Since the Mars 
Observer PMIRR experiment will obtain temperature profiles ex- 
tending from the surface up to about 80 km altitude, it will be 
extensively sampling middle atmospheric regions in which gravity 
wave drag may play a dominant role. Estimating the drag then 
becomes crucial to the estimation of the atmospheric winds from the 
PMIRR-observed temperatures. 

An iterative diagnostic model based upon one previously devel- 
oped and tested with Earth satellite temperature data [4] will be 
applied to the PMIRR measurements to produce estimates of the 
small-scale zonal drag and three-dimensional wind fields in the 
Mars middle atmosphere. This model is based on the primitive 
equations, andean allow for time dependence (the time tendencies 
used may be based upon those computed in a Fast Fourier Mapping 
procedure). The small-scale zonal drag is estimated as the residual 
in the zonal momentum equation; the horizontal winds having first 
been estimated from the meridional momentum equation and the 
continuity equation. The scheme estimates the vertical motions 
from the thermodynamic equation, and thus needs estimates of the 
diabatic heating based upon the observed temperatures. The latter 
will be generated using a radiative model. It is hoped that the 
diagnostic scheme will be able to produce good estimates of the 
zonal gravity wave drag in the Mars middle atmosphere, estimates 
that can then be used in other diagnostic or assimilation efforts, as 
well as in more theoretical studies. 

References: [1] Barnes J. R. (1990) JGR, 95, 1401-1421. 
[2] Theodore B. et al. (1993) JAS, submitted. [3] Theodore B. et al. 
(1993 ) Icarus, submitted. [4] Marks C. J. (1989) 7A£ 46, 2485-2508. 


EDDY MIXING COEFFICIENT VALUES DERIVED FROM 
SIMULATIONS WITH THE NASA AMES MARS GCM, 
J. R. Bamcs and T. Walsh, College of Oceanic and Atmospheric 
Sciences, Oregon State University, Corvallis OR 97331 -2209, USA. 

Values of eddy mixing coefficients, especially that for vertical 
mixing, are of particular importance for various photochemical 
models of the Mars atmosphere, which are either globally orzonally 
averaged. These models represent atmospheric transport processes 
in terms of eddy diffusion. While this is not appropriate for the 
adveedve transport by the mean meridional circulation (in general, 
a part of the eddy diffusion tensor actually does represent advec- 
tion — a ''correction” to the Eulerian-mean circulation), it can be 
applied to simulate both the small- and larger-scale eddy transports 
of an atmospheric constituent 

An effort is underway to estimate values of the eddy mixing 
coefficients for the Mars atmosphere using circulation data gener- 
ated with the NASA Ames Mars GCM. This model simulates the 
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three-dimensional winds in the atmosphere, and these are then used 
as inputs for a tracer transport model. The latter model has previ- 
ously been used in interactive dust transport experiments with the 
Ames GCM [1], Idealized tracer transport experiments, with a 
conservative or nearly conservative hypothetical tracer, are per- 
formed to generate data from which the eddy diffusion coefficents 
can be estimated. Carrying outmatched pairs of tracer experiments, 
using two very different initial states, permits all four components 
of the diffusion tensor to be determined [2]. A large number of GCM 
experiments have been conducted, spanning a range of seasons and 
atmospheric dust loading, allowing the eddy mixing coefficients to 
be estimated for a variety of atmospheric conditions. 

References: [1] Murphy J. R. (1993) JGR, 98, 3197-3220. 
[2] Plumb R. A. and Mahlman J. D. (1993) JAS, 44, 298-327. 


MARS DUST-DRIVEN TIDES AND THEIR IMPACT ON 
THE THERMOSPHERE. S. W. Bougheri and R. W. Zurek 2 , 
1 Lunar and Planetary Laboratory, University of Arizona, Tucson AZ 
85721, USA, 2 Jet Propulsion Laboratory, Pasadena CA 91109, 
USA. 

It has been known since the early Mariner 6, 7, and 9 missions 
that dust loading of the lower atmosphere, and the subsequent 
aerosol heating during dusty periods, impacts the martian middle 
and upper atmospheres. A quantitative measure of this lower atmo- 
sphere forcing was obtained by the Viking 1 and 2 landers, from 
which observed amplitudes of semidiurnal surface-pressure oscilla- 
tions were correlated with normal-incidence dust optical depths 
[1—3]. It appears that the dominant semidiurnal mode is a good 
indicator of global dust content or mean dust optical depth, espe- 
cially during dust storm events. A classical tidal model that repro- 
duces the surface pressure oscillations measured by these Viking 
landers in 1977 was used to calculate tidal amplitudes and phases 
up to -43 km [3]. These tidal characteristics were calculated for 
various dust optical depth conditions ranging from typical dusty 
periods to global dust storm times. Reasonable extrapolations can 
be made to higher altitudes if one assumes that the vertically 
propagating tidal modes continue to grow without dissipation or 
breaking. It is very likely that gravity waves also play an important 
role in the structure and dynamics of the middle atmosphere of Mars, 
since the large topographical relief should produce vigorous gravity 
wave fluxes [4], 

Semidiurnal tidalmodes, significantly enhanced by lower atmos- 
phere dust-induced heating, may indeed propagate to the Mars 
thermosphere (£100 km) before breaking and generating turbu- 
lence. The preferential enhancement of the semidiurnal tides during 
dust storm onset is primarily due to the elevation of the tidal heating 
source in a very dusty atmosphere. The (2,2) semidiurnal tidal mode 
was shown to have the largest variation with dust optical depth, as 
measured by Viking landerinstruments [3]. Also, the (2,2) mode has 
the largest vertical wavelength of all the semidiurnal tidal modes, 
and thus is most likely to penetrate into the thermosphere before 
breaking and to modify the largely in situ solar-driven behavior 
otherwise expected. The tides may also be partially responsible for 
determining the height of the martian homopause (-125 km). 

Dust-driven tides and their vertical propagation into the lower 
thermosphere are also important because tidal effects may explain 
the unexpected behavior of the Mars thermospheric temperatures 


. and densities observed during the Mariner 9 nominal and extended 
' missions [5,6]. This was a period when Man experienced a raging 
global dust storm and its gradual abatement Solar fluxes were also 
near a minimum in the 11 -year periodic cycle. The Mariner 9- 
inferred exospheric temperatures of -325 ± 40 K are significantly 
warmer than the 250 K value given by the solar-driven simulation 
of a three-dimensional thermospheric general circulation model 
(MTGCM) using EUV-UV forcing alone [6,7]. This is consistent 
with an expanded (warmer) lower and middle atmosphere that is 
reflected in an ionospheric peak altitude nearly 20-30 km higher 
than expected. A further discrepancy between model and observed 
fields is evident in the Mariner 9-infcrred O -fractional abundances 
as a function of local time. The calculated O variation is much 
different than inferred from UVS 1 30.4-nm airglow intensity scans 
across the dayside. This discrepancy indicates an additional forcing 
of the dynamics that is responsible for redistributing O about the 
planet [5,6]. Evidently, the 1971 dust storm had a significant impact 
upon the Mariner 9-observed temperatures and densities. These 
outstanding features imply a significant lower-upper atmosphere 
coupling that modifies the otherwise diurnal Mars thermospheric 
structure driven in situ by the Sun. . 

TheNadonalCenterfor AtmosphericResearch (NCAR) MTGCM 
code is the most recent in a hierarchy of thermospheric models 
developed to examine the thermospheric structure and dynamics of 
Earth, Venus, and Mars. The MTGCM code is based on the pres- 
sure-coordinate primitive equations of lower atmosphere dynamic 
meteorology. The physical processes incorporated into this model, 
however, are those appropriate to thermospheric dynamics, i.e., fast 
molecular vertical diffusion of heat, momentum, and constituents at 
thermospheric heights . The MTGCM model presently sol ves primi- 
tive equations yielding global solutions for the zonal, meridional, 
and vertical velocities, total temperatures, geopotential heights, and 
important neutral (C0 2 , CO, O, N 2 , NOX) and dayside ion (OJ, 
COJ, 0+, NO+) densities on a 5° latitude and longitude grid. The 
vertical coordinate is log -pressure, with 24 levels spanning roughly 
100-250 km. A vertical spacing of two grid points per scale height 
is achieved. The boundary conditions are zero fluxes at the top and 
prescribed winds, densities, temperatures, and the geopotential at 
the bottom. 

A prototype MTGCM tidal model was recently developed to 
incorporate the effects of the upward propagating semidiurnal tides 
by modifying the MTGCM lower boundary conditions near 100 km 
according to classical tidal theory [6,8]. The basic variable used to 
parameterize the lower boundary is the geopotential. It is formulated 
as a function of latitude (via Hough functions), longitude, and time 
in units of days. The other model fields are then derived from the 
geopotential using the standard classical tidal relations [8]. The 
lower boundary is thus specified by the complex coefficients, in- 
cluding the amplitude and phase, of the geopotential for each mode. 
Specific lower boundary tidal amplitudes and phases for the (2,2) 
semidiurnal mode were estimated from Zurek’s classical tidal 
model [2,3]. Values were calculated for a moderate value of dust 
optical depth (t = 1.5) corresponding to average dusty conditions 
near southern summer solstice. This prescription effectively couples 
the lower and upper atmospheres in a physically realistic manner. 

MTGCM simulations incorporating these semidiurnal tides show 
a dramatic change in the divergence and convergence of horizontal 
and vertical winds, thereby modifying the diurnal temperature and 
density structure throughout the Mars thermosphere. The inclusion 
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of the (2,2) semidiurnal mode in the model gives rise to enhanced 
midaftemoon temperatures that are closer to the values inferred 
from the Mariner 9 UVS (see Figs. 1 and 3). The diurnal variation 
of atomic O is also shown to be strongly dependent on the specified 
tidal amplitudes and phases, i.e., a greatly improved match of 
Mariner 9 UVS and modeled O mixing ratios is obtained (see 
Fig. 2). Increased attenuation of the semidiurnal propagating tides 
exists at high solar activity (similar to Earth), rendering the upper 



Fig. 1. Solar cycle variation of MTGCM exospheric temperatures as a function 
of local tune about the equator. Mariner 6, 7 and 9 solar only (a) and tidal cases 
(b) are illustrated. All figures taken from [6]. 



Fig. 2. Diurnal variation of the atomic O mixing ratio at the height of the 
dayside ionospheric peak. Mariner 6. 7, and 9 solar only (a) and tidal cases (b) 
are illustrated. The Mariner 9 tidal model distribution more closely matches that 
observed by the Mariner 9 UVS. 


thermospheric response weakened with respect to that at solar 
minimum. Equivalently, tidal effects are discemable higher in the 
thermosphere during solar minimum, since the dissipation due to 
eddy/molecular viscosity and thermal conductivity is diminished 
(see Fig. 3). This suggests that dust-driven tides, especially during 
solar minimum time periods, cannot be ignored when addressing the 
Mars thermospheric structure and dynamics. 

Further study of the coupling of the Mars lower and upper 
atmospheric regions requires the development of a tidal model that 
explicidy addresses the assumptions of tidal amplitude growth and 
regular phase shifting throughout the middle atmosphere (50- 
100 km). The development of a single coupled chemical-energy- 
dynamical model that realistically addresses the important physics 
in the entire martian atmosphere is extremely difficult. Various 
timescales for chemical, radiative, and dynamical processes change 
radically with altitude. However, it is feasible to run two well- 
developed general circulation models in tandem: the upper bound- 
aiy of one would be used to specify the lower boundary of the other. 
Work is underway to modify the NASA Ames GCM model [e.g., 
9,10] to cover 0-80 km; the NCAR MTGCM code will be extended 
to a 70-km lower boundary. The overlap will permit the Ames GCM 
to be run for various lower-atmosphere dust conditions and seasons. 
The GCM pressure surface near 70 km will then be used to progres- 
sively extract tidal mode amplitudes and phases for specification of 
the MTGCM geopotential height field. This procedure will provide 
an improved way to examine the propagation of dust-storm-driven 
tides and their effects throughout the entire atmosphere from the 
ground to the exobase. These two models will provide an improved 
set of tools for the interpretation of future PMIRR temperature data 
from the Mars Observer mission. 



Fig. 3. Mariner 9 MTGCM tidal simulation: alatitude slice of total temperature 
at22.5°S, with contours from 1 10-280 K. Midafternoon exospherictemperarures 
are 280 K, much warmer than without the (2,2) tidal mode. 
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CLOUD, DUST, AND OZONE VERTICAL PROFILES 
FROM SOLAR OCCULTATION MEASUREMENTS: 
IMPLICATIONS FOR DYNAMICS. E. Cbassefiere and J. E. 
Blamont, Service d' Aeronomie du CNRS, BP 3, 91371 Venieres- 
le-Buisson Cedex, France. 

An instrument has been designed for solar occultation measure- 
ments of the martian atmosphere from the Phobos spacecraft. It was 
composed of three different dispersive systems working in the 
ultraviolet (UV: 0.22-0.32 pm) for the measurement of 0 3 and 
aerosols, in the near infrared (NIR: 0.76 pm, 0.94 pm) for the 
detection of 0 2 and H 2 0, and in the infrared (JR: 1.9 pm, 3.7 pm) 
where C0 2 and H 2 0 were measured. A detailed description of the 
instrument may be found in the special issue of Nature [1], Its 
principle objective is to measure from the Phobos orbit the spectrum 
of the Sun, modified by atmospheric extinction, during sunset. The 
UV-NIR spot has an angular diameter of 1 aremin, or »3-km vertical 
resolution, and is located near the center of the solar disk. The 1R 
field is about twice as large and its line of sight is shifted by 
8 5 aremin, or *20 km, relative to the previous one. It is therefore 
located near the edge of the solar disk. S ampling times are generally 
0.5, 1, and 2 s for IR, UV, and NIR channels respectively, corre- 
sponding to vertical excursions of the line of sight of 1, 2, and4km 
respectively under nominal conditions. The instrument operated 
from February 8 to March 26 (the martian equinox occurred on 
February 17). The latitude of the intersection of the Sun-spacecraft 
axis with the surface of Mars varied from -1 1° to +20°, the seasonal 
date L, being in the range 0°-20°. All measurements were therefore 
made near northern spring equinox in equatorial regions. Due to an 
error in the pointing system, only partial results were obtained, the 
region below *30 km altitude being never sounded by the UV-NIR 
spectrometer. On the contrary, nine complete occultations were 
obtained in the 1R channels, whose line of sight was fortunately 
*20 km below the UV-NIR axis. 

Five clouds were detected above 45 km altitude and their vertical 
structure recorded at six wavelengths between 0.28 pm and 3.7 pm 
[2]. They have a small vertical extent (3-6 km) and a vertical optical 
depth less than 0.03. The thermal structure, as derived from satu- 
rated profiles of water vapor observed by the instrument in the 
infrared [3] , does not allow the C0 2 fros t point to be reached at cloud 
altitude, strongly suggesting that cloud particles are formed of H 2 0 
ice. Under the assumption of spherical particles, a precise determi- 
nation of their effective radius, which varies from cloud to cloud and 
with altitude, was obtained, ranging from 0.15 to 0.85 pm. An 
estimate of the effective variance of the particle size distribution is 
*0.25 . The number density of cloud particles at the peak extinction 
level is *1 cm* 3 . Assuming that particles are levitated by eddy 
mixing, the eddy diffusion coefficient, K, was found to be 10 5 — 


lO 6 cm 2 s* 1 at 50 km by comparing the results of a simple steady- 
state model with observations. Indeed, since the clouds areobserved 
at twilight, which coincides with the diurnal maximum of the 
ambient temperature, they can be assumed to be in a steady state. If 
their thermodynamic state were to vary quickly during the day, our 
optical thickness at twilight would correspond to unrealistic values 
in earlier hours when the temperature is lower. An atmospheric 
temperature of 165-170 K at *50 km is inferred. The negative 
temperature gradient above the cloud is large (1.5-2 Kkm- 1 )* 

Dust was also observed and monitored at two wavelengths, 1 .9 
and 3.7 pm, on nine different occasions [2]. The top of the dust 
opaque layer, defined as the level above which the atmosphere 
becomes nearly transparent at the wavelengths of observation, is 
located near 25 km altitude, with variations smaller than ±3 km from 
place to place.The scale height of dust at this altitude is 3-4 km. The 
effective radius of dust particles near the top of the opaque layer is 
0.95 ± 025 pm and increases below with a vertical gradient of 
*0.05 pm km- 1 . Assuming that particles are levitated by eddy 
mixing, the eddy diffusion coefficient, K, is found to be *10® cm 2 s- 1 
at 20-25 km. Ah effective variance of 025 (±50%) for the dust size 
distribution is obtained on the basis of a simple theoretical model for 
the observed vertical gradient of the effective radius of dust 

Concerning ozone, only six occultations may be analyzed in the 
near-UV wing of the Hartley continuum (*0.270-0.283 pm). Ozone 
is detected at one time with a high degree of confidence below 
*60 km and above *37 km altitude [4]. Its number density shows a 
maximum of *10* cm ~ 3 around 42—45 km. Its scale height above 
45 km is 6 ± 2 km. Ozone, present in two other cases with a similar 
vertical profile, is not detected, thatis, notin excess of «2x 10 7 cm- 3 , 
for the three other observations. An analysis based on a one- 
dimensional steady-state photochemical model indicates that when 
ozone is present the eddy diffusion coefficient, K, is *2 x 10* cm 2 
s^ 1 around 40-45 km altitude within a factor of *2 on either side. 
Below 50 km, the water vapor mixing ratio is lower than *2 x 10~ 5 . 
Such a low ratio is consistent with a temperature decreasing with 
increasing altitude, with a typical temperature of *150 K at *50 km 
altitude. Above 50 km, the H 2 0 mixing ratio is definitely lower than 
10- 5 , strongly suggesting saturation at high altitude. The agreement 
between the modeled and observed profiles of ozone is significantly 
improved by assuming a sharp increase of K over a typical scale of 
1 atmospheric scale height around 40—45 km altitude, from 
*5 x 1G 5 cm 2 s' 1 below 40 km up to «10 7 cm 2 s* 1 above 45 km. The 
absence of ozone found in at least half the cases may be simply 
reproduced by using the thermal structure and humidity conditions 
derived from the analysis of the clouds [2], i.e., a temperature and 
a H 2 0 mixing ratio of respectively *165 K and “10- 4 at 50 km 
altitude. 

The first implication of these measurements is the good agree- 
ment between the different estimates of the eddy diffusion coeffi- 
cient K. From a simple comparison between observed number 
density and size profiles of water ice aerosols and dust particles, the 
value of K is nearly the same at low and high altitude, i.e., 
*10* cm 2 s- 1 around 20-25 km and lOMOScm 2 s- 1 around 50 km. 
From ozone measurements, a most likely value of 2 x 10 6 cm 2 s -1 
within a factor of 2 on either side is found around 40-45 km, 
consistent with values derived from cloud and dust measurements. 
The suspected sharp increase of Kby more than 1 order of magnitude 
between 40 km and 45 km around the pivot value of 2 x 1 0 6 cm 2 s* 1 
could be the signature of breaking gravity waves. From a two- 
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dimensional modeling of the effect of breaking gravity waves on the 
general circulation and temperature field in the middle atmosphere 
of Mars [5] using Lindzen's parameterization, which explains well 
millim eter wave observations made with the IRAM 30-m telescope 
[6], such an increase, by about 1 order of magnitude per half 
atmospheric scale height, may be easily reproduced around 50 km 
altitude, in the breaking region. Because the ozone profile is ex- 
tremely sensitive to a sharp increase of K around 30-50 km, the 
future observation of ozone by stellar and solar occultation from the 
Mars 94 spacecraft (SPICAM experiment) could bring some inter- 
esting information on the breaking level of internal gravity waves. 

The second implication concerns the likely presence of an inver- 
sion of the temperature profile around 35 Ion when clouds arc 
observed, in order to allow the rather high temperature (165-170 K) 
and thermal gradient (-1.5-2 K km* 1 ) inferred at the cloud level. 
More precisely, the presence of two condensation levels at *15 km 
and *50 km altitude are inferred from water vapor and cloud 
measurements respectively. Above *15 km, a regular decrease of 
the water vapor mixing ratio up to *35 km is observed [3] and it is 
interpreted as the signature of a saturated atmosphere. Around *50 
km, thin clouds, very likely formed of water ice, are detected [2]. 
This means that above *35 km exists a region where the atmosphere 
is not saturated in H 2 0, containing enough water vapor to give rise 
to the formation of clouds. The reason for the temperature inversion 
might be adiabatic wanning due to convective motions induced by 
breaking gravity waves. It was shown from the two-dimensional 
modeling that weak temperature inversions are generated above the 
equatorial region at solstice for certain parameters of the waves [5], 
and a similar phenomenon could be produced by the model at 
equinox if the preliminary results of the last millimeter- wave obser- 
vations made at IRAM arc confirmed (Lellouch, personal commu- 
nication). The second issue is how water molecules are transported 
from the low atmosphere up to the cloud region. A possible expla- 
nation is that upwelling in the ascending branch of the Hadley cell 
is able to supply water ice in the shape of small particles, across the 
cold saturated region (*15-35 km), to the nonsaturated layer below 
the clouds, where ice may sublimate and condense again at cloud 
level (*50 km). A simple heuristic model of water ice vertical 
transport at equator has been constructed [7]. Different constraints 
have been introduced relating to the coagulation probability during 
the upward transport phase, the minimum amount of H 2 0 that must 
be deposited at the sublimation level, and the infrared optical 
thickness of ice particles. It was shown that these constraints are 
jointly satisfied for reasonable values of parameters, indicating that 
a quantity (zonally averaged) of *0.0 1/ pm of water ice in the shape 
of small particles of radius 0.1-0.2 pm might be present in the 
middle atmosphere (15-35 km) above equatorial regions. 

The third implication is that, from a qualitative point of view, 
thin clouds and an ozone layer seem to be the signatures of two well- 
defined states of the middle atmosphere, the region around40-50 km 
being warm and wet when clouds are observed (15% of cases), cold 
and dry when ozone is present (*15-50% of cases). An outstanding 
fact is that, considering these two typical atmospheric states, the 
vertical structure of the observed constituent, i.e., water ice particles 
and ozone molecules respectively, as well as the inferred atmo- 
spheric parameters (temperature, water vapor mixing ratio, eddy 
diffusion coefficient), do not seem to change significantly from 
place to place, or from time to time. It means that the variability of 
the middle atmosphere of Mars could not consist in a wide and 


continuous range of possible states, with random fluctuations of 
temperature and water vapor amount, but rather in an organized 
system of a small number of stable configurations. This general 
result reinforces our belief that the thermal structure and the degree 
of humidity of the middle atmosphere of Mars are driven by dynami- 
cal processes. The thermal structure could be controlledby meridional 
winds induced by breaking gravity waves, with associated inversion 
of the meridional and vertical thermal gradients. Eddy mixing 
partially reflects the interaction between gravity waves and general 
circulation, and a sharp change of K at some altitude, as suspected 
from our analysis of ozone, could be the local signature of this 
interaction. It must be noted that an intermediate state, for which 
there is an abundance of water vapor too small to yield the formation 
of clouds but too large to prevent the ozone density from reaching our 
detection threshold around 40 km, is also observed (*35-70% of 
cases), and the previous schematic picture will certainly have to be 
refined. The simple typology established from Phobos 2 solar occul- 
tation measurements (see Table 1 ) could be the starting point toward 
a more extensive study of the dynamical processes controlling the 
middle atmosphere of Mars with the help of oncoming observations 
by the Mars Observer mission (PMIRR, Radio Science experiment). 


TABLE 1. 


Water vapor mixing ratio 

<2xl(>5 

2-6x10-5 

6x10-5 

Type of atmosphere 

cold/dry 


warm/wet 

Ozone 

yes 

no 

no 

Cloud 

no 

DO 

yes 

Percentage of occurrences 

(15-50%) 

(35-70%) 

(15%) 
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PARTICLE SIZES AND COMPOSITION OF MARS 
ATMOSPHERIC DUST BASED UPON VIKING AND MARI- 
NER 9 OBSERVATIONS. R. T. Clancy 1 , S. W. Lee 1 , and G. R. 
Gladstone 2 , 'Laboratory for Atmospheric and Space Physics, 
University of Colorado, Boulder CO 80309, USA, 2 Southwcst 
Research Institute, P.O. Drawer 285 10, 6220 Culebra, San Antonio 
TX 78228, USA. 

Mars atmospheric dust can play an important role in the thermal 
structure of the Mars atmosphere during periods of high dust loading 
[1]. However, the radiative properties of Mars atmospheric dust 
remain uncertain due to uncertain definitions of the dust composi- 
tion and size distribution. The analysis by Toon et al. [2] of Mariner 
9 IRIS spectra during the 1971-1972 global dust storm indicated a 
reasonable match between the modeled 9-pm absorption of 
montmorillinite and the observed 9-pm absorption. Toon et al. also 
determined that an effective (cross-section weighted) mean radius 
of 2.5 pm (Rniode = 0*4 Mm) provided a consistent fit of montmorillinite 
to the IRIS dust spectra at 9 pm. Pollack et al. [1] analyzed Viking 
lander observations of atmospheric extinction and scattering at 
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visible-near IR wavelengths (0.5-1 .Opm), and obtained consisten- Method: The model solves the primitive equations on the 

cy with the Toon et al. dust size distribution when the effects of ) ' sphere using a spectral method in the horizontal and finite differ- 
nonspherical particle shapes were included. An additional, minor ences in the vertical. The coarse resolution (-5° longitude x 
(1%) component of visible-ultraviolet absorbing material was re- 5° latitude and 10 vertical levels) requires some parameterization of 
quired to model the derived visible (0.86) and ultraviolet (0.4-0.6) subgridscale friction and heating. It is tempting to embaik on a 
[3] single-scattering albedos of the dust, since montmorillinite does scries of complicated representations of radiation, turbulence, etc. 
not absorb sufficiently in this wavelength region [4], This soon becomes a burden with the specification of alarge number 

There are, however, several deficiencies with the Toon et al. of arbitrary parameters. Also the computational cost increases, 
model of Mars atmospheric dust As Toon et al. and subsequent which precludes long integrations and regime studies. We have 
analyses of Viking IRTM observations [5] pointout, montmorillinite developed two distinct models: ( 1) The Simplified Global Circula- 
exhibits strong, structured absorption near 20 pm that is not exhib- tion Model (SGCM) has for its thermodynamic forcing a linear 
ited by the IR observations. Furthermore, the observed ratio of relaxation to a prescribed radiative equilibrium temperature profile 
visible-to-infrared extinction opacities for Mars atmospheric dust is on a certain timescale. Frictional forces are similarly linear and are 
approximately 2 [6], which is twice the calculated ratio for mont- applied at the lower boundary of the model. (2) The Intermediate 
morillinite with a cross-section weighted mean radius of 25 pm Global Circulation Model (IGCM) has a prescribed surface temper- 
[7,8]. In addition, higher values for the visible single-scattering ature distribution and a boundary layer scheme that mixes heat and 
albedo of the dust (0.92 vs. 0.86) have been derived from the Viking momentum into the atmosphere. Radiative cooling is assumed 
IRTM albedo channel, as well as a corrected value of 0.55 for the constant throughout the atmosphere and a dry convective adjust- 
single-scattering asymmetry parameter [8]. ment scheme is in place. 

W c have conducted a combined analysis of the Viking IRTM and Both models provide a computationally cheap and conceptually 

Mariner 9 observations to reassess the model of Mars atmospheric simple framework in which to study regimes of planetary scale flow, 

dust in terms of a single optical component, which might satisfy the No Mountain Flows: The Viking landers relayed data that 

ultraviolet- to-mfrared measurements of dust absorption and scatter- suggested the existence of baroclinic waves in the northern winter 
ing. The optical constants for palagonite [9,10] are incorporated in midlatitudes. Among the seemingly random fluctuations there ap- 
a doubling- adding radiative transfer model of the Mars atmosphere peared highly periodic phases when synoptic systems passed over 
to simulate Mariner 9 IRIS spectra as well as the Viking IRTM IR the landers with monotonous regularity [1]. This would be in stark 
band observations. We also derive visible and ultraviolet single- contrast to the Earth’s rather chaotic weather patterns. Martian 
scattering albedos based on the Hansen and Travis Mie scattering baroclinic waves may be closer to the regular waves found in the 
code [11]. A tentative conclusion is that smaller dust particles rotating annulus experiment [3], 

(Rmode = 0.15 pm, cross-section weighted mean R = 1.2 pm) com- In order to study such baroclinic waves the SGCM was set up to 

posed of palagonite provide a much improved fit to the Mariner 9 run with no mountains and hence remove the possibility of station- 

IRIS spectra; agreement with the observed ratio of visible-to-infra- ary waves. Radiative equilibrium temperatures, representing south- 

red extinction opacities; and ultraviolet and visible single-scatter- em hemisphere winter, were taken from [5]. The sensitivity of the 
ing albedos comparable to their observed values. modeled climatology was then investigated by varying (principally) 

References: [1] Pollack J.B.etal.(1979)/GR, 84, 2929— 2945. the frictional and thermal relaxation timescales. 

[2] Toon O. B. et al. (1977) Icarus , 30, 663-696. [3] Pang K. et al. In general the modeled baroc lini c waves appear as small undu- 
(1976) Icarus, 27, 55—67. [4] Pollack J. B. et al. (1977) JGR, 82, lations propagating on a zonally orientated jet stream (Fig. la). The 
4479-4496. [5] Hunt G. E. (1979) JGR, 84, 8301-8310. [6] Martin fast radiative relaxation timescale of the atmosphere limits the 
T. Z. (1986) Icarus, 66, 2—21. [7] Zurek R. W. (1982) Icarus, 50, amplitude of the waves. Compare this with the Earth’s midlatitudes 
288-310. [8] Clancy R. T. and Lee S. W. (1991) Icarus, 93, where cyclones grow to finite amplitude in turbulent storm tracks. 
135—158. [9] Roush T. (1993), personal communication. [10] Clark In many of the integrations, however, a zonal wavenumber 3 pattern 
R. N. et al. (1990) JGR, 95, 14463-14480. [11] Hansen J. E. and (Fig. lb) emerges and grows to large amplitude, dominating the 
Travis L. D. (1974) Space Sci Rev,, 16, 527-610. flow. This wave 3 propagates eastward around the midlatitudes at 

a constant phase speed for at least 1000 days (the longest run). The 
resulting flow is very reminiscent of the single wavenumber steady 

ASPECTS OF THE MARTIAN GLOBAL CIRCULATION, waves found in the rotating annulus laboratory experiments [3]. 

M. Collins andI.N. James, Department of Meteorology, University Read [6] and White [7] have suggested that the steady annulus 

of Reading, Reading RG2 6AU, UK. waves may be “free modes” of the full nonlinear governing equa- 

tions. It is assumed that, in the interior of the fluid, heating and 
Introduction: Recent advances in high-speed computers have friction are small and the nonlinear terms are zero, i.c., if u.V£ is a 

given atmospheric scientists a new laboratory in which to work, nonlinear term u. = 0 while neither u nor V£ is zero. A similar 
Numerical models of the Earth's atmosphere and its dynamics have approach is used in modon studies of, e.g., atmospheric blocks and 
provided great insight. Given the lack of observations of the martian jovian vortices. The assumption of zero heating is clearly not valid 
atmosphere it is natural to apply such models to help us understand here. Also the existence of the wave 3 is sensitive to both heating and 
aspects of, for instance, large-scale circulation patterns. Here we frictional timescales. A close examination of the potential vorticity 
have modified a global circulation model of the Earth's atmosphere budgetreveals that (away from the ground) a stable balance is struck 
for use with martian parameters and have studied the relevant between the nonlinear and heating terms, the rapid radiative timescale 
regimes of planetary flow. allowing the wave to be maintained at finite amplitude. 
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Integrations with other radiative equilibrium temperature pro- 
files have large amplitude zonal wavenumber-4 patterns (Fig. 2) 
with similar properties, so we believe these results to be robust. The 
connection between these waves and the regular oscillations found 
in the Viking data is unclear, especially considering the relative 
simplicity of the SGCM. It is desirable then to add an extra degree 
of realism to the model. 

Flows with Mountains: The SGCM uses a terrain-following 
C coordinate system in the vertical. In these simulations the more 
recent DTM orography is used. Again the southern hemisphere 
winter radiative equilibrium temperatures from Pollack et al. [5] 
were interpolated onto the model grid and the flow was spun up from 
a resting isothermal state. 

Orography generates stationary waves that can cither propagate 
vertically through the atmosphere or exponentially decay away from 
the surface. Vertically propagating waves can mix heat andmomcn- 
tum meridionally and thus contribute to the maintenance of the 
global circulation. Partitioning eddy statistics into stationary and 
transient components showed that the stationary waves were domi- 
nating with the transients virtually absent. A similar result has been 
found by [2]. This would constitute a considerable asymmetry 
between the hemispheres, the northern winter midlatitudes being 
the only place where traveling synoptic systems are present 

The eddy kinetic energy maximum for these integrations is near 
the model top (Fig. 3), which is a rigid lid. This implies that 
stationary waves arc being spuriously reflected when really they 
should propagate out of the model domain. To partially resolve this 
problem linear friction (sponge layer) is added to “soak up” the 
energy. Examination of the eddy statistics then showed the exist- 
ence of some transient behavior. It appeared that the stationary 
waves were saturating and choking off the transients. In other 
integrations, however, the sponge layer reduced the transient ampli- 
tude. This shows that a careful treatment of the top boundary of any 
model of the martian atmosphere is necessaiy. 
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Fig. 2. 


The IGCM: Although the SGCM provides adequate simula- 
tions of the martian circulation, some aspects of the flow are not well 
represented. The introduction of more complicated physical 
parameterizations requires the specification of more arbitrary con- 
stants. These are not certain even for the Earth and are difficult to 
determine for Mars. Here, we have kept the level of parameteriza- 
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Fig. 4. 


tion to a minimum. A boundary layer scheme with prescribed 
surface temperatures heats the atmosphere, which is in turn cooled 
at a constant rate. Convective instability is possible and a simple 
adjustment scheme is applied. The model reproduces fairly well the 
zonally averaged climatology of [2] (Fig. 4). 

Examining the flow at the Viking lander sites in the model shows 
the existence of mobile baroclinic systems. A spectral analysis of 
pressure has strong spectral peaks at similar frequencies to those 
found by [1]. The flow shows high temporal coherence compared to 
that of the Earth. 

An interesting aspect that was First pointed out to us by [4] 
concerns the annual pressure cycle. The large annual signal in the 
Viking data has been attributed to the melting and freezing of the 
polar C0 2 frost caps and would constitute a significant change in the 
mass of the atmosphere. This model has no evaporation/condensa- 
tion and the atmospheric mass is conserved. A sizable annual 
pressure cycle is, however, still evident. The redistribution of mass, 


required by geostrophic balance of the midlatitude jets, gives a large 
pressure gradient term that fluctuates throughout the year as the 
flow switches from easterly to westerly. 

Conclusions: The two models described have succeeded in 
modeling some observed aspects of the martian atmospheric circu- 
lation and some aspects that have been generated by other, more 
complicated GCMs. The regularity of synoptic systems on Mars sets 
it aside from the Earth where weather systems are highly unpredict- 
able. Martian baroclinic waves may be more akin to the regular 
waves found in rotating annuli. The large amplitude orography also 
gives a different regime with stationary waves swamping the flow. 
This is in contrast to the Earth’s northern hemisphere where station- 
ary and transient waves are of equal importance. 

References: [1] Barnes J. R. (1981) J. Atmos . Sci., 38, 225- 
234. [2] Barnes J. R. et al. (1993) JGR, in press. [3] Hide R. and 
Mason P. J. (1975) Adv. Phys ., 24, 47-100. [4] Hourdin F. (1992) 
Ph.D. thesis, I’UniversitS Paris 7. [5] Pollack J. B. et al. (1981) J. 
Atmos . Sci., 38, 3-29. [6] Read P. L. (1985) Dyn. Atmos. Oceans, 
9, 135-270. [7] White A. A. (1986) Q.J.R. Met. Soc ., 112, 749-774. 


IMPACT OF THE C0 2 AND H 2 0 CLOUDS OF THE 
MARTIAN POLAR HOOD ON THE POLAR ENERGY 
BALANCE. F. Forget and J. B. Pollack, NASA Ames Research 
Center, Moffett Field CA 94035-1000, USA. 

Clouds covering extensive areas above the martian polar caps 
have regularly been observed during the fall and winter seasons of 
each hemisphere. These “polar hoods” are thought to be made of 
H 2 0 and C0 2 . In particular, the very cold temperatures observed 
during the polar night by Viking and Mariner 9 around both poles 
have been identified as C0 2 clouds [1,2] and several models, includ- 
ing GCM [3], have indicated that the C0 2 can condense in the 
atmosphere at all polar latitudes. 

Estimating the impact of the polar hood clouds on the energy 
balance of the polar regions is necessary to model the C0 2 cycle and 
address puzzling problems like the polar caps assymetry. For ex- 
ample, by altering the thermal radiation emitted to space, C0 2 
clouds alter the amount of C0 2 that condenses during the fall and 
winter season. 

We have analyzed the complete set of Viking IRTM data to 
define the spatial and temporal properties of the polar hoods, and 
estimated how their presences affect the energy radiated by the 
atmosphere/caps system to space. 

The IRTM observations provide good spatial and temporal cov- 
erage of both polar regions during fall, winter, and spring, when a 
combination of the first and the second Viking year is used. Only the 
IRTM brightness temperatures at 1 1, 15, and 20 pm are reliable at 
martian polar temperatures. To recover the integrated thermal fluxes 
from the IRTM data, we have developed a simple model of the polar 
hood, thought to consist of “warm” H 2 0 clouds lying above colder 
and opaque C0 2 clouds. Such a model is based on the analysis of the 
IRIS spectra [4], and is consistent with the IRTM data used here. 

The high thermal emissivity of the C0 2 ice particles is possible 
because of contamination by dust and water [5], Besides, we have 
found new evidence for such opaque C0 2 clouds in the IRTM data. 

The model has been calibrated and validated by testing its ability 
to reproduce polar hood IRIS spectra from simulated IRTM tem- 
peratures. 
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A specific treatment of the IRTM data has also allowed us to 
recover the “true” ground-ice temperatures, which only depend on 
the atmospheric pressure (solid-gas equilibrium), but which were 
difficult to measure directly because of the presence of the C0 2 and 
H 2 0 clouds. 

We have also obtained estimations of the amount of water ice in 
the polar clouds during winter, and a climatology of the frequency 
and temperature of the C0 2 clouds. 

Thus, we have been able to estimate the wavelength-integrated 
thermal emission by the surface/atmosphere system at the top of the 
atmosphere and to compare this with what would have been emitted 
by the bare polar caps. In the future, this information will allow us 
to address polar process questions with a realistic parametrization 
of the polar hood. 

References: [1] HuntG. E. (1980) GRL, 7, 481-484. [2] Paige 
D. A. (1985) Ph.D. thesis, Caltech. [3] Pollack J. B. et al (1990) 
JGR, 95, 1447-1473. [4] Paige D. A. (1992) personal communica- 
tion. [5] Warren S. G. et al. (1990) JGR, 95, 717-741. 


MULTIANNUAL SIMULATIONS WITH THE MARS 
CLIMATE MODEL- R. M. Haberle 1 , H. C. Houben 2 , and R. E. 
Young 1 , ! NAS A Ames Research Center, Moffett Field CA 94035- 
1000, USA, 2 Space Physics Research Institute, Sunnyvale CA 
94087-1315, USA. 

We have simulated the martian atmospheric circulation for 10 
Mars years using a simplifed three-dimensional climate model that 
is briefly described below. Analysis of these simulations reveals 
several interesting behaviors not previously recognized in other 
three-dimensional models. The first concerns the potential for 
interannual variability in the martian climate system. We now know 
from the aperiodic frequency of planet-encircling storms, as well as 
differences in atmospheric water vapor from one year to the next, 
that the real planet exhibits interannual variations. The results of 
our simulations also show significant interannual variability and, 
interestingly, that variability maximizes during southern spring, the 
dust storm season. A second and related finding is that transient 
eddies in the southern hemisphere during winter arc indeed weaker 
than their northern hemisphere counterparts as was found in more 
sophisticated models. However, southern hemisphere eddy activity 
rises dramatically during spring and is greater than that in the 
northern hemisphere at any season. 

The model we use is based on the primitive equations in log- 
pressure coordinates. The model is spectral (spherical harmonics) 
in the horizontal, and finite-differenccd in the vertical. For the 
multiannual simulations, wc run with trapezoidal truncation using 
6 zonal waves, and a maximum total wavenumber of 16. There are 
14 vertical levels that gradually increase in thickness from -500 m 
in the lowest layer to about one-half scale height in the upper layer. 
The model runs on a dedicated Silicon Graphics Workstation; with 
half-hour time steps, it takes about 8 hours to simulate 1 Mars year. 

However, the difference between this model and a GCM is its 
greatly simplified physics. We use Newtonian cooling with a 2-day 
damping time for thermal forcing, with a Rayleigh friction “sponge” 
layer at the top and a quadratic drag force at the bottom. The rest of 
the atmosphere is inviscid. The relaxation temperature field is 
zonally symmetric and is computed offline from a separate model. 


In spite of these simplifications the model zonal mean and eddy 
statics are very similar to those of more sophisticated models. 

To assess the nature of the model’s interannual variability we 
have looked at the hemispherically integrated zonal (ZKE) and eddy 
(EKE) kinetic energies per unit mass. In the southern hemisphere 
EKE is lowest during summer (-10 mks units) and then steadily 
rises during fall. However, just before winteritpeaks (-35) and then 
declines to a local minimum (-25) at solstice. It remains near this 
minimum until mid winter, when it rises, dramatically reaching 
-100 by spring. Following this there is an abrupt decline in EKE, 
which appears to be related to a weakening in the wave-3 distur- 
bances. The decline is short-lived, however, as EKE rises dramat- 
ically once more but this time to even higher values (-125) by mid 
spring. This is the largest value of EKE in either hemisphere at any 
time of year. The resurgence in EKE at this time is not well 
understood, but it occurs at a time when ZKE is rapidly declining 
and wave-1 disturbances arc rapidly growing. This suggests that 
barotropic processes are playing a role, but more analysis is needed 
to confirm this. 

Interestingly, the EKE in the northern hemisphere shows a 
similar seasonal pattern, although the amplitude of the fluctuations 
is much less pronounced. Furthermore, the springtime peak in 
northern hemisphere EKE (-50) is comparable to that during late 
fall and mid winter, which is considerably different from that in the 
southern hemisphere. 

We have also looked at the maximum and minimum EKE for any 
given seasonal date for the 10-ycar simulation. This we take to be 
a measure of the magnitude of interannual variability. The greatest 
variability occurs in the southern hemisphere during the late winter 
and spring seasons. At this time, EKE can be as high as 150 and as 
low as 50. Thus, southern-hemisphere EKE can vary from year to 
year by as much as a factor of 3. Evidently, this variability is also 
manifest in the model’s “surface stress,” which also shows similar 
variations. 

We find these results to be very intriguing, particularly since they 
occur in a model with simple but repeatable forcing. They suggest 
that the martian climate system has enough nonlinearity to produce 
substantial interannual variability. 

REGARDING TRACER TRANSPORT IN MARS’ WINTER 
ATMOSPHERE IN THE PRESENCE OF NEARLY 
STATIONARY, FORCED PLANETARY WAVES. J. L. 
Hollingsworth 1 , R. M. Haberle 1 , and H. C. Houben 2 , J NASA Ames 
Research Center, Moffett Field CA 94035, USA, 2 Space Physics 
Research Institute, Sunnyvale CA 94087, USA. . 

Large-scale transport of volatiles and condensates on Mars, as 
well as atmospheric dust, is ultimately driven by the planet’s global- 
scale atmospheric circulation. This circulation arises in part from 
the so-called mean meridional (Hadley) circulation that is associ- 
ated with rising/poleward motion in low latitudes and sinking/ 
equatorward motion in middle and high latitudes [1]. Intimately 
connected to the mean circulation is an eddy-driven component due 
to large-scale wave activity in the planet’s atmosphere. During 
winter this wave activity arises both from traveling weather systems 
(i.c., barotropic and baroclinic disturbances) [2] and from “forced” 
disturbances (e.g., the thermal tides and surface-forced planetary 
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waves) [3,4]. In this paper we investigate possible contributions to 
the effective (net) transport circulation from forced planetary waves. 

A three-dimensional global spectral primitive equations model 
is applied in this investigation to simulate “perpetual” northern and 
southern winter solstice conditions. Previously this model has been 
applied to investigate the sensitivity of winter transient-eddy activ- 
ity to changes in thermal dissipation and to examine the transport 
circulation associated with traveling baroclinic eddies [5,6]. Re- 
cently the model has been modified to include kinematic effects of 
topography using a terrain-following vertical coordinate (i.e., effec- 
tively log-o, where a - p/p s ). In the calculations to be presented, a 
trapezoidal truncation of total wavenumber 16 and zonal wave- 
number 6 is used and simplified physics are applied: Thermal 
radiative forcing is prescribed in the form of constant Newtonian 
cooling, and momentum drag is specified in the form of a height- 
dependent Rayleigh friction. Planetary-scale surface zonal asym- 
metries in Mars’ topography excite planetary waves having deep 
vertical structure. Because of the nature of the mean circulation and 
the planet’s size, forced responses can be “global” (occurring in both 
eastern and western hemispheres) and can become “focused” into a 
finite latitude band in the winter extratropics, similar in respects to 
results of a linear analysis of forced stationary waves on Mars [4]. 
In the presence of a nonzero mean meridional circulation, we evalu- 
ate the effective steady transport circulation associated with the 
planetary-wave activity and compare it with the so-called “residual” 
mean meridional circulation, which is frequently used as an ap- 
proximation to the former. In particular, we examine potential 
hemispheric asymmetries in the net transport circulation due to 
differences in the wave forcing in northern and southern hemi- 
spheres. We also investigate changes to the effective transport under 
highly dusty, northern winter conditions. Potential impacts steady 
planetary waves can have on the poleward transport of tracers 
during northern and southern winter will be discussed. 

References: [1] Haberle R. M. et al. (1993) JGR , 98 , 
3093-3124. [2] Barnes J. R. et al. (1993) JGR , 98, 3125-3148. 
[3] Zurek R. W. and Haberle R. M. (1988) 7. Atmos . Sci., 45, 
2469-2485. [4] Hollingsworth J. L. and Barnes J. R. (1993) 
J . Atmos. Sci., submitted. [5] Haberle R. M. et al. (1993) in prep- 
aration. [6] Houben H. et al. (1989) Eos, 70, 1179-1180. 


ANALYZING MARTIAN WINDS AND TRACER CON- 
CENTRATIONS USING MARS OBSERVER DATA, 
H. Houben, Space Physics Research Institute, Sunnyvale CA 94087- 
1315, USA. 

During the course of a day, the Mars Observer spacecraft will 
acquire globally distributed profiles of the martian atmosphere. It is 
highly desirable that this data be assembled into synoptic weather 
maps (complete specifications of the atmospheric pressure, tem- 
perature, and winds at a given time), which can in turn be used as 
starting points in the study of many meteorological phenomena. 
Unfortunately, the special nature of the Mars Observer data pre- 
sents several challenges above and beyond the usual difficult prob- 
lem of data initialization. 

Mars Observer atmospheric data will consist almost exclusively 
of asynoptic vertical profiles of temperatures (or radiances) and 
pressures, whereas winds are generally in balance with horizontal 


gradients of these quantities (which will not be observed). It will 
therefore be necessary to resort to dynamical models to analyze the 
wind fields. As a rule, data assimilation into atmospheric models 
can result in the generation of spurious gravity waves, so special 
steps must be taken to suppress these. In addition, the asynoptic 
nature of the data will require a four-dimensional (space and time) 
data assimilation scheme. The problem is to find a full set of 
meteorological fields (winds and temperatures) such that, when 
marched forward in time in the model, they achieve a best fit (in the 
weighted least-squares sense) to the data. 

The proposed solution is to develop a model especially for the 
Mars Observer data assimilation problem. Gravity waves are fil- 
tered from the model by eliminating all divergence terms from the 
prognostic divergence equation. This leaves a diagnostic gradient 
wind relation between the rotational wind and the temperature field. 
The divergent wind is diagnosed as the wind required to maintain 
the gradient wind balance in the presence of the diabatic heating. 
The primitive equations of atmospheric dynamics (with three prin- 
cipal dependent variables) are thus reduced to a simpler system with 
a single prognostic equation for temperature — the variable that will 
be best observed. (This balance system was apparently first derived 
by Chamey as a first-order Rossby number expansion of the equa- 
tions of motion.) Experience with a full primitive equation model of 
the martian atmosphere indicates that a further simplification is 
possible: at least for short-term integrations, the model can be 
linearized about the zonally symmetric basic state. 

The above linearized spectral model has the property that the 
solution at any point (in space and time) is a known linear function 
(represented by a matrix operator) of the initial state. Finding the 
initial state that best fits the observed data requires an inversion of 
this matrix. If the problem is overdetermined (more data points than 
unknowns), the best fit in the least-squares sense is obtained by a 
generalized inverse (singular value decomposition) of the matrix. 
The proposed model has 14 vertical levels, 7 zonal waves, and 16 
meridional modes, for about 1500 complex variables. The PMIRR 
instrument alone, with approximately 1000 vertical profiles per 
day, should easily provide enough data in a day to overdetermine the 
model. Thus, a daily weather map appears to be quite feasible. (Of 
course, it is probable that there are errors in this model. There is, 
however, a generalization of the singular value decomposition — 
the total least-squares solution — that allows for errors in the model 
as well as in the observations.) There are two classes of motion in 
the martian atmosphere that are not well represented by the gradient 
wind relation — atmospheric tides and equatorial Kelvin waves. 
These modes must be independently modeled and removed from the 
observations before the assimilation procedure is performed. (The 
fact that Mars Observer will orbit at a nearly constant time of day 
will inhibit the direct observation of atmospheric tides in any case.) 

Atmospheric tracers, such as water vapor and dust, are easily 
incorporated into the balance system model with continuity equa- 
tions. The sole difficulty arises from the possibility that a spectral 
tracer model will show negative tracer amounts at various points. It 
is therefore necessary to constrain the solution of the least-squares 
problem so that the tracer concentrations at all of the model’s (-620) 
transform grid points are positive. This is facilitated by the introduc- 
tion of surface source/sink terms into the model. Again, there arc 
mathematical algorithms available for the solution of this problem. 
The above program indicates that the assimilation of Mars Observer 



1 6 Workshop on Atmospheric Transport on Mars 


profiles into daily synoptic maps of temperature, winds, and tracers 
is quite feasible, if computationally demanding. 

THE ANNUAL PRESSURE CYCLE ON MARS: RESULTS 
FROM THE LMD MARTIAN ATMOSPHERIC GENERAL 
CIRCULATION MODEL. F. Hourdin, F. Forget, and 
O. Talagrand, Laboratoire de Meteorologie Dynamique du CNRS, 
Ecole Normale Superieure, Paris, France. 

At Laboratoire dc Meteorologie Dynamique (LMD) in France 
we have been developing a General Circulation Model (GCM) of 
the martian atmosphere since 1989. The model has been described 
rather extensively elsewhere [1] and only the main characteristics 
are given here. 

The dynamical part of the model, adapted from the LMD terres- 
trial climate model, is based on a finite -difference formulation of the 
classical “primitive equations of meteorology.” The radiative trans- 
fer code includes absorption and emission by C0 2 (carefully validat- 
ed by comparison to line-by-line calculations [2]) and dust in the 
thermal range, and absorption and scattering by dust in the visible 
range. Other physical parameterizations are included: modeling of 
vertical turbulent mixing, dry convective adjustment (in order to 
prevent vertical unstable temperature profiles), and a multilayer 
model of the thermal conduction in the soiL Finally, the condensa- 
tion-sublimation of C0 2 is introduced through specification of a 
pressure-dependent condensation temperature. The atmospheric 
and surface temperatures are prevented from falling below this 
critical temperature by condensation and direct precipitation onto 
the surface of atmospheric C0 2 . The only prespecified spatial fields 
are the surface thermal inertia, albedo, and topography. 

The first simulations were performed for a visible dust optical 
depth constant in both space and time, x Va = 0.2, a value typical of 
clear-sky conditions [3]. Starting from an isothermal state at rest, 
without polar caps, the model quickly adjusts (in less than one 
martian year) to a regime in which the annual atmospheric mass 
cycle, due to condensation of the atmospheric C0 2 into the polar 
caps, is essentially repeatable from one year to the next This mass 
cycle is in great part responsible for the large-amplitude seasonal 
oscillations of the surface pressure recorded by the Viking landers. 
Those oscillations were rather well reproduced in those first simu- 
lations, as well as the synoptic pressure oscillations observed during 
winter that arc due to traveling planetary waves [1]. The agreement 
with temperatures retrieved from the Mariner 9 infrared measure- 
ments is also very good. 

However, the surface pressure cycle was found to vary strongly, 
with latitude. This meteorological contribution, due to intemalmass 
redistribution, can be divided into two parts: the first one results 
from the modulation by temperature of the latitudinal pressure 
variations due to orography; the second one is a purely dynamical 
contribution due to geostrophic balance between the strong zonal 
winds and latitudinal pressure gradients [1]. These dynamical and 
orographic contributions appear to have comparable magnitudes 
and are stronger in high latitudes where the dynamical effect be- 
comes predominant At 75°S, for instance, the dynamical compo- 
nent contributes to a 25% temporal variation of the surface pressure 
over the course of the year, about the same magnitude as the 
variations due to the condensation-sublimation cycle. The two 
effects tend to reinforce each other in midlatitude and counteract 


each other in polar regions. Longitudinal variations also strongly 
amplify the orographic effect 

More recent simulations were performed in order to test the 
sensitivity of the model to the most uncertain parameters and 
especially to the visible albedo of the ice, A^, and its thermal 
emissivity, e^. Model results have shown that they strongly influ- 
ence the energy and mass atmospheric budget [3,4]. In our simula- 
tions, the best agreement to Viking pressure data is obtained for very 
different values of these parameters for the two caps. Using A ice = 
0.62 and = 0.6 for the northern cap and A ice = 0.45 and e icc = 0.7 6 

for the southern cap, very good agreement is achieved as shown in 
Fig. 1 . The latitudinal time evolution of the edge of both caps is also 
very close to Viking orbiter data as shown in Fig. 2 (except at the 
very end of the recession phase). As already mentioned by other 
authors [4] these “best fit” values of the ice parameters must not be 
interpreted as real physical properties of ice materials since they 
also account for some model deficiencies and, in particular, for the 
nonrepresentation of C0 2 ice clouds, resulting from atmospheric 



Fig. 1. A comparison of the time evolution of the pressure simulated at both 
Viking landing sites to that observed (*1?est-fit” simulation). The thin curves are 
the observed pressure for year 1 (dotted), 2 (solid) and 3 (dashed) of the mission. 
The black heavy curve is the simulated pressure, which has been smoothed by 
applying a 15-sols average. Sols are martian solar days; sol 0 corresponds to the 
landing of Viking l. 



sols 


Fig. 2. The time evolution of the mean latitude of the polar cap boundaries as 
simulated by the LMD GCM and deduced, during cap recession, from Viking 
orbital imaging. 
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condensation during the cap formation. This may explain the very The amount of dust activity in the northern hemisphere is clearly 

low value of the "best fit” emissivity of the northern cap, which has very relevant to the role of atmospheric transport in the dust cycle, 
been observed to be preferentially covered by polar haze. The low The classic global storms that occur during spring in the southern 
emissivity will account for the trapping of thermal radiation by hemisphere are observed to transport dust from sources in the 
clouds during the formation phase. However, this value of emissiv- southern hemisphere to sinks or temporary depositories in the 
ity is too low for the recession phase during which no atmospheric north. The question of whether atmospheric transport can close the 
condensation occurs and would produce too fast a recession of the dust cycle, Le„ return the dust to the southern hemisphere sources 
cap. This is compensated by a high value of the “best fit” albedo on some timescale, is clearly relevant to the solution of the p uzzl e 
(which has a larger impact during the recession than during the of how the dust storm cycle is modulated, i.e., why storms occur in 
formation phase since condensation essentially occurs in the polar some years but not in others. 

night), which decreases the recession rate by reducing the part of the There are data that suggest that the spring/early summer season 

incoming solar radiation absorbed by the ice. in the northern hemisphere of Mars during the year following the 

It must be noted that the dust amount has been kept at its clear- major 1977 storms observed by Viking was very dusty. A number 

sky value Tyj, = 0.2, whereas atmospheric dust content was ob- of observations of the vicinity of the receding north polar cap 

served to strongly increase during northern winter, even during the showed clear evidence of substantial dust activity in the sub- Arctic 

clear years of the Viking mission. This strong simplification could region. 

be relaxed in the future, either by specifying the dust atmospheric Although some of the activity was less well organized and less 
content as a function of time and space (for instance by using the extensive than what was observed near the south cap during spring, 
results of the IRTM experiment aboard Viking [5] or results from at least one set of observations define a regional dust storm event 
fiituremartian missions), or by modeling the atmospheric transport, that, if merely an example of dust storms in the northern hemi- 

which has not been done in the LMD GCM until now. sphere, could play a crucial role in closing the dust cycle. Though 

As far as the atmosphere is concerned, an important scientific by- the observations have been previously reported [2], their possible 
product of future missions would be an accurate estimation of the significance relative to the question of atmospheric dust transport 
atmospheric mass budget The evaluation of the meteorological justifies a review in this workshop. Observations on three consecu- 
contribution to the temporal and spatial pressure variations will also five orbits of Viking Orbiter 2 at Ls = 65 arc consistent with a single 
be very useful since it contains indirect but accurate information on storm moving to the south -southeast (prograde and equatorward) at 
the near-surface winds and may help to validate the circulation an average speed of 6 m s- 1 . The storm is observed in the region to 
models. Finally, the spatial results of future martian missions will the north of Olympus Mons and Alba Patera between 60° and 75° 
be fully exploited only if they arc used in conjunction with an latitude and between 150° and 110° longitude, 
adequate description of the physical laws of the atmosphere, as Unfortunately there was no color coverage of this event, so no 
provided by an atmospheric general circulation model. This is definitive conclusions on composition can be reached. The designa- 

cuirentlydoneforweatherforecastingusingdataassimilationmeth- tion of this storm as a dust storm is based on its very low contrast 
ods. For Mars, since we are concerned with a posteriori analysis, the with the surface and with the morphology, which is similar to other 

problem of the computational cost of the assimilation techniques is dust storms observed on the planet. The dust storm is not mono- 

much less crucial than for operational forecasting and some power- lithic; rather, the storm is composed of many discrete billows 
ful methods such as variational assimilation [6] may easily be used, similar in structure to those observed within the 1977a global storm 

The development of those methods for Mars is presently under study observed by Viking. 

at LMD [7]. The leading edge of the storm is defined by a very sharp front. 

References: [1] Hourdin F. et al. (1993) J. Atmos . Sci. r in and the trailing edge was also well defined. Since the coverage on 

press. [2] Hourdin F. (1992) JGR, 97, 18319-18335. [3] WoodS. E. any single orbit of the spacecraft was insufficient to image the entire 
and Paige D. A. (1992) Icarus, 99, 1—14. [4] Pollack J. B. ct al. leading or trailing edges of the storm, only a lower limit on the 
(1993) JGR, 98, 3149-3141. [5] Martin T. Z. (1986) Icarus, 66, dimension transverse to the motion can be made; the largest ob- 
2—21 . [6] Talagrand O. (1992) In Energy and Weather Cycles in the served dimension is roughly 500 km, but the geographical distribu- 
Chmate System (A. Raschke and D. Jacob, eds.), 187-213, NATO tion of the three images suggests a considerably larger extent The 
Advanced Study Institute. [7] Talagrand O. and Hourdin F. (1990) dust is more narrowly confined to a strip 350km wide parallel to the 
International Symposium on Assimilation of Observations in Me - direction of motion with a strong suggestion of expansion between 
teorology and Oceanography, 578—580, World Meteorological the earliest and latest views. 

Organization. With only one example in the class, the frequency of this sort of 

northern regional dust storm is clearly undetermined Similarly, the 
question of whether these are annual events or whether they take 

NORTHERN HEMISPHERE DUST STORMS ON MARS, place only in years following global storms is completely open. It is 
P. B. James, University of Toledo, Toledo OH 43606, USA. hoped that the synoptic monitoring provided by the MOC experi- 

ment will be able to answer the first question, though the answer 
Dust storms in the northern hemisphere of Mars appear to be less will still be confined to only a single year with its own unique dust 

common than the more familiar southern hemisphere storms, and history. 

essentially no activity north of about 30° latitude has been docu- References: [1] Martin L J. and ZurekR. W. (1993) JGR, 98, 
mented [1]. The data are, however, subject to an observational bias 3221-3247. [2] James P. B. (1985) Recent Advances in Planetary 

because Man is near aphelion during oppositions, which occur Meteorology, Cambridge, 85-100. 
during the most likely seasons for dust activity in the north. 


18 Workshop on Atmospheric Transport on Mars 


THE INFLUENCE OF OROGRAPHY ON THE TRANS- 
PORT OF ATMOSPHERIC CONSTITUENTS. M. M. Joshi, 
S. R. Lewis, and P. L. Read, Atmospheric, Oceanic and Planetary 
Physics, Clarendon Laboratory, Parks Road, Oxford 0X1 3PU, UK. 

Given the fact that martian orography varies over a range of two 
scale heights (cf., about one scale height on Earth), it is thought that 
on Mars some aspects of the circulation should be affected to a large 
degree by the presence of orography. Initial results from two numeri- 
cal models are presented that illustrate the effect of orography on 
cross-equatorial transport on Mars. 

On Earth, western boundary currents are found not only in the 
oceans as one would expect, but also in the troposphere, the most 
notable example being the East African Jet, which has an important 
effect on the monsoon circulation of the Eastern Indian Ocean [1]. 

Two models were used for this investigation: a simple GCM 
(hereafter referred to as the SGCM) and a simple one-level baro tropic 
model (BM). The SGCM, originally developed by Hoskins and 
Simmons [2], solves the hydrostatic primitive equations on a sphere 
by representing model variables by spherical harmonics in the 
horizontal, with nonlinear terms being calculated in grid-point 
space. In the vertical, the model uses the a coordinate system in 
finite-difference form. The model represents the radiative processes 
in the atmosphere by a Newtonian relaxation to a specified zonal 
mean temperature state, and surface drag by Rayleigh friction in the 
lowest model level. 

BM is a nonlinear spectral barotropic model [3] , which solves the 
barotropic vorticity equation on the sphere 

D£/Dt = -f ■ V • V (h/H) + (£o- QH (1) 

Co represents a specified zonal mean absolute vorticity state to 
which the model is relaxed over a timescale T. h is the height of the 
orography and H is a scale height. The main advantage of this model 
is that many runs can be carried out relatively quickly and runs can 
also be carried out at spatial resolutions, which would be impractical 
with the SGCM. 

Nonlinear barotropic boundary current theory states that, in the 
steady state, the width of an inertial (i.e., where friction is negli- 
gible) boundary current will be given by (U/p) 1/2 . It is to be noted, 
however, that while the width of the jet is given by (U/p) 1/2 , the 
width of the core of the jet, where most of the mass is transported, 
is much less than this [1] and so higher model resolution, as well as 
significandy changing the character of the orography that controls 
the jet, may actually be needed to resolve the core of the jet. 

As an initial experiment, BM is initialized with a restoration 
corresponding to a constant small zonal velocity representative of 
that near the surface. The steady-state velocity vectors, overlaid on 
an orography contour plot, are shown in Fig. 1 . The run is carried out 
at resolution T21, equivalent to a grid-point spacing of about 6°. It 
can be seen that there is a tendency for the cross-equatorial flow to 
form into a western boundary current, especially at about 90°W, i.e., 
at the eastern edge of the Tharsis plateau. The width of this current 
is approximately (U/P) 1/2 . The effect of model resolution on this 
result can be tested by repeating this experiment at T42 or even T63. 



Fig. 1. Longitude/latitude plot or velocity vectors superimposed on orography 
contours (interval 1km): (U/b) 1/2 scale is shown at bottom, right. 
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Fig. 2. Comparisons of cross-equatorial mass transport for SGCM model runs 
at T21 and T42. 


While BM can show regions where the effect of orography is 
likely to be important, the SGCM is needed to show directly how 
orography affects mass transport as it is only in the SGCM that the 
mountains act explicitly as boundaries and not just as large vorticity 
forcing regions. Initial results compare favorably with the BM and 
do tend to show the presence of western boundary currents. 

Initial SGCM results also suggest that runs at model resolution 
T21 give slightly less intense western boundary jets when compared 
to runs done atT42, as shown in Fig. 2. T21 runs, therefore, may not 
sufficiently resolve these boundary currents and hence not ade- 
quately model cross-equatorial mass transport 

The presence of these boundary currents in the tropics will have 
important implications for mass transport in the lower branch of the 
Hadley Cell on Mars, and, as this cell plays such an important role 
in the martian circulation, on the cross-equatorial transport of 
atmospheric constituents such as dust and water vapor, it is impor- 
tant to model them as accurately as possible. 

References: [1] Anderson D. L. T. (1980) Orographically 

Controlled Cross-Equatorial Flow, G.A.R.P. Publications Series, 
23, 317-355. [2] Hoskins B. J. and Simmons A. J. (1975) Q. J. R. 
Met. Soc., 101, 637-655. [3] James I. N. (1988) Q. J. R. Met. Soc., 
114, 619-637. 




LPl Technical Report 93-05 19 


TEMPORAL AND SPATIAL MAPPING. OF SURFACE 
ALBEDO AND ATMOSPHERIC DUSf' OPACITY ON 
MARS. S. W. Lee 1 , R. T. Clancy 1 , and G. R. Gladstone 2 , 
^boratoiy for Atmospheric and Space Physics, University of 
Colorado, Boulder CO 80309, USA, 2 Southwest Research Institute, 
P.O. Drawer 28510, 6220 Culebra, San Antonio TX 78228, USA. 

The Mariner 9 and Viking missions provided abundant evidence 
that eolian processes are active over much of the surface of Mars 
[1,2]. Past studies have demonstrated that variations in regional 
albedo and wind-streak patterns are indicative of sediment transport 
through a region [3,4], while thermal inertia data [derived from the 
Viking Infrared Thermal Mapper (IRTM) dataset] are indicative of 
the degree of surface mantling by dust deposits [5-9]. The visual and 
thermal data are therefore diagnostic of whether net erosion or 
deposition of dust-storm fallout is taking place currently and whether 
such processes have been active in a region over the long term. These 
previous investigations, however, have not attempted to correct for 
the effects of atmospheric dust loading on observations of the 
martian surface, so quantitative studies of current sediment trans- 
port rates have included large errors due to uncertainty in the 
magnitude of this “atmospheric component” of the observations. 

We have developed a radiative transfer model that allows the 
atmospheric dust opacity to be determined from IRTM thermal 
observations. Corrections for the effects of atmospheric dust loading 
on observations of surface albedo can also be modeled. This ap- 
proach to determining “dust-corrected surface albedo” incorporates 
the atmospheric dust opacity, the single-scattering albedo and par- 
ticle phase function of atmospheric dust, and the bidirectional 
reflectance of the surface, and accounts for variable lighting and 
viewing geometry. The most recent dust particle properties [10,1 1] 
are used. The spatial and temporal variability of atmospheric opac- 
ity (t) strongly influences the radiative transfer modeling results. 
This approach allows t to be determined at the highest spatial and 
temporal resolution supported by the IRTM mapping data; maps of 
“dust-corrected surface albedo” and atmospheric opacity can be 
constructed at a variety of times for selected regions. Information on 
the spatial and temporal variability of surface albedo and atmos- 
pheric opacity, and inferences of the amount of dust deposition/ 
erosion related to such variability, result. 

Analyses of IRTM mapping observations of the Syrtis Major 
region, covering a time span of more than a martian year, will be 
presented. 

References: [1] Veverka J. etal. (\971)JGR, 82, 4167-4187. 
[2] Thomas P. etal. (1981) Icarus, 45, 124-153. [3] Lee S. W. et al. 
(1982 )JGR, 87, 10025—10042. [4] Lee S. W. (19S6) Symposium on 
Mars: Evolution of its Climate and Atmosphere (V. Baker et al., 
eds.), LPI Tech. Rpt. 87-01, 71-72. [5] Kieffer H. H. et al. (1977) 
JGR, 82, 4249-4295. [6] Christensen P. R. (1982) JGR, 87, 9985- 
9998. [7] Christensen P. R. (1986) JGR, 91, 3533-3545. 
[8] Christensen P. R.(1986)/carwj, 68, 217-238. [9] Jakosky B. M. 
(1986) Icarus, 66, 1 1 7- 1 24. [ 1 0] Clancy R. T. and Lee S. W. ( 1 99 1 ) 
Icarus , 93, 135-158. [1 1] Clancy R. T. et al., this volume. 


DYNAMICS OF THE ATMOSPHERE OF MARS: REVIEW 
OF SOME OUTSTANDING PROBLEMS. C. B. Leovy, 
Department of Atmospheric Sciences, University of Washington, 
Seattle WA 98195, USA. 

Theoretical and modeling studies based on fragmentary data are 
responsible for our current level of understanding of the dynamics 
of Mars' atmosphere at the planetary scale. Largely because the 
boundary conditions and governing parameters are such that the 
large-scale circulation regime is similar to that of Earth, it has been 
possible to derive a remarkable amount of information from a 
relatively small number of observations [1]. Nevertheless, the level 
of detail available from theoretical and modeling studies far exceeds 
the validation capability of existing data. Global observations from 
Mars Observer instruments are expected to go far to remedy this 
mismatch. This review will present and examine several outstand- 
ing dynamical questions, which Mars Observer measurements should 
help to answer. These include: What are the temporal and three- 
dimensional spatial distributions of thermal tides? What is the 
distribution of momentum flux divergence due to internal gravity 
waves? Are normal modes an important component of the general 
circulation? If so, how do they vary with season? Do normal modes 
contribute to the initiation of planetaiy-scale dust storms? What are 
the details of the global circulation response to large and rapid 
injections of dust? Current understanding of these problems, the 
potential role of Mars Observer observations for addressing them, 
and the implications of answers for our understanding of the present 
climate and past climate regimes will be discussed. 

References: [1] Zurek etal. (1992) in Mars (H. Kieffer et al., 

eds.), 835-933, Univ. of Arizona. 


MODELING AND DATA ASSIMILATION FOR MARS 
OBSERVER. S. R. Lewis and P. L. Read, Atmospheric, Oceanic 
and Planetary Physics, Department of Physics, University of Oxford, 
Parks Road, Oxford OX1 3PU, UK. 

The Mars Observer mission will return a substantial set of 
meteorological data for the atmosphere of Mars, principally in the 
form of thermal infrared soundings from the Pressure Modulator 
Infrared Radiometer (PMIRR). The dataset will enable global stud- 
ies of a range of dynamical and physical phenomena in the martian 
atmosphere. The polar orbital configuration of Mars Observer, 
however, means that the coverage of the martian atmosphere and 
surface by remote sounding instruments will be asynchronous. This 
leads to difficulties of interpretation when studying synoptic-scale 
phenomena with timescales of a few days. In an attempt to overcome 
such problems we propose to use data assimilation techniques, as 
currently employed for operational weather forecasting on the Earth, 
in conjunction with a Martian General Circulation Model (MGCM) 
that is under development at Oxford and Reading Universities. An 
overview of the model and assimilation scheme will be given, with 
reference to problems encountered so far in applying such methods 
to remote-sensing observations of the martian atmosphere. 

The MGCM is based on a semi spectral primitive equations 
model [1], which covers the surface to 80 km altitude, using about 
20 o levels in the vertical, and with a triangular spectral truncation 
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at wavenumber 42 in the horizontal, equivalent to a grid spacing 
of 2.8° in latitude and longitude. The resolution is chosen to be 
somewhat greater than that necessary to resolve the smallest fea- 
tures potentially visible in the PMIRR data in order to model 
realistically the effects of topography, etc. The martian surface is 
represented using the U.S. Geological Survey Digital Topographic 
Map (DIM). S ome experiments with a simple form of the model are 
described in a companion submission (Joshi et al.). Parameteriza- 
tion schemes for poorly understood processes arc kept simple wher- 
ever possible. New schemes for theMGCM are underdevelopment 
in collaboration with a small group at Reading University. These 
include a boundary layer scheme, a dry convection scheme, a diffu- 
sive soil model, and a scheme for orographic gravity wave drag. The 
radiation scheme will initially be based on that developed at 
Laboratoire de Meteorologie Dynamique, Paris [2], with the possi- 
bility of further refinement to include the effects of nonlocal thermo- 
dynamic equilibrium in the upper atmosphere in collaboration with 
colleagues at Granada University. 

Data assimilation will be conducted using a modified form of the 
Analysis Correction Scheme that is in current operational use at the 
U.K. Meteorological Office [3]. This scheme is an approximation to 
the optimum analysis, which involves a form of successive-correc- 
tion iterated between each time step of the model. Continuous 
(asynoptic) repeated insertion is conducted sequentially, assimilat- 
ing each measurement at the time and place it is valid without prior 
interpolation. Observational increments to the model are spread in 
time and space by empirically determined functions with scales that 
depend on the difference between the current model time and the 
time at which the data are valid. Observations typically enter the 
model five hours before they become valid with a low weight and 
with their influence spread over large horizontal scales. As an 
observation’s valid time approaches, its weight increases and the 
scale it affects becomes smaller. Finally, an observation is discarded 
about one hour after the model has passed its valid time. The weight 
of each observation is scaled by the relative magnitudes of observa- 
tional and forecast model error. Within each cycle of the assimila- 
tion scheme each model variable is analyzed sequentially followed 
by the derivation of multivariate increment fields for dynamical 
balance in order to avoid losing much of the newly introduced 
information. For example, following the analysis of the temperature 
field from PMIRR soundings, the thermal wind equation will be 
used to derive velocity increments in balance with the temperature 
increments (Fig. 1). 

The Analysis Correction Scheme has been found to produce 
improved analyses in the Southern Hemisphere on Earth compared 
to the previous Optimum Interpolation scheme in use at the U.K. 
Meteorological Office, primarily because observational increments 
arc spread over larger areas. In the case of Mars, however, there will 
be no surface observations and the use of satellite temperature 
soundings alone may introduce particular problems. Surface pres- 
sure measurements are strictly required for direct information about 
the barotropic component of the circulation, and for verification of 
the mass of the model atmosphere, which otherwise will depend on 
the highly uncertain parameterizations of atmospheric C0 2 conden- 
sation and sublimation in polar regions. Estimates of surface pres- 
sure may be retrieved from PMIRR nadir observations and will be 
used by the assimilation scheme whenever possible. Radio occulta- 



Flg. L Assimilation scheme cycle with multivariate balancing increments. 

tion measurements may also be used to derive more accurate surface 
pressures at a few locations per sol at high latitudes; these measure- 
ments will also be assimilated if available. 

Trials of the scheme are currently being conducted using syn- 
thetic data generated from parallel model runs. In future it will be 
important to test the assimilation scheme with artificial data gener- 
ated by a different model in order to avoid the “identical twin” 
problem. Once the Mars Observer spacecraft is in its mapping orbit 
and verified retrievals of PMIRR data are routinely available (around 
early 1994) the assimilation scheme and model may be run in batch 
mode, processing about one week’s data in a day on a dedicated 
workstation. Typical output from the assimilation will consist of 
fields as represented in the model (velocity, temperature, and sur- 
face pressure) stored at regular intervals from which derived quan- 
tities and diagnostics (potential vorticity, energy budgets, transport 
of heat, momentum and material tracers, Eliassen-Palm fluxes, etc.) 
may be calculated and which may be used, for example, as initial 
states for further independent model studies. 

Major scientific objectives to be addressed during the course of 
the mission will include (1) determination of the role of transient 
baroclinic waves in heat and momentum transport on Mars, (2) the 
factors influencing their spatial and temporal regularity, and (3) the 
influence of surface orography on mass transport and of other major 
meteorological variables and tracers (see also the companion con- 
tribution by Joshi et al.). It is also hoped to make use of the zonally 
averaged assimilated transport fields in conjunction with retrieved 
fields of water vapor from PMIRR in studies of the synoptic and 
seasonal transport of water during the Mars Observer period. 

References: [1] Hoskins B. J. and Simmons A. J. (1975) Q. J. 
R . Met . Soc., 101, 637-655. [2] Hourdin F. (1992) JGR, 97, 
18319-18335. [3] Lorenc A. C. et al. (1991) Q. J . R. Met . Soc ., 117, 
59-89. 
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MARS OPERATIONAL ENVIRONMENTAL SATELLITE To understand the transport of tracers in the martian atmosphere, 

(MOES)— A POST-MARS OBSERVER DISCOVERY * and particularly to identify their sources and sinks, it is necessary to 
MISSION. S. S. Limaye, University of Wisconsin-Madison, have systematic global, synoptic observations that have yet to be 
Madison WI 53706, USA, and the MOES Team. attained. Clearly these requirements are not easy to achieve from a 

single spacecraft in orbit, but if we focus on specific regions of the 
The study of the general circulation of planetary atmospheres in planet, e.g., polar vs. low and mid latitudes, then it is possible to 
terms of the constituent atmospheric transport processes has been a attain a nearly ideal coverage at a reasonable spatial and temporal 
subject just beyond our reach for most planets, despite the numerous resolution with a system of just two gatwlljt^ Mars Observer is 
spacecraft missions. The primary reason has been the lack of ad- about to yield good coverage of the polar latitudes, so we focus 
equate systematic spatial and temporal observations. We do not yet initially on the region not covered well in tw™. of diurnal coverage, 
have any systematic observations of importantmartian atmospheric and in terms of desired observations. These observations will pro- 
phenomena such as regional and global dust storm evolution and vide the initial data for the numerical models of the martian weather 
decay stages, planetary-scale waves, and cloud systems. Yet such andclimatethatcanbeverifiedonly with better temporal and spatial 
systems are key components of the processes invoked to explain data. 


major features in the evolution of the martian environment such as 
the polar layered terrain and regional albedo changes. Mars Ob- 
server (MO) presents the first real opportunity of making global 
atmospheric and surface observations, but its Sun-synchronous 
polar orbit limits the temporal and global coverage and resolution. 
Observations ±60° latitudes are limited to wi thin about an hour of 
2:00 a.m. and 2:00 pm. local time throughout the mission life. Thus 
while the diurnal tidal amplitudes will be well determined (supple- 
mented by modeling), the precise phase relationships of the diurnal 
components and the amplitudes andphascs of the semidiurnal tides 
and other higher-frequency waves and the associated astmospehric 
transports will not be sampled. At polar latitudes the high-inclina- 
tion orbit provides substantial overlap so that the polar processes 
will indeed be well observed in local time over the mission duration, 
but their connection to the lower latitudes will remain largely 
unknown. 

Lacking an infrared imaging system, MO can sample only a 
fraction of the martian globe with nighttime observations of the 
atmospheric phenomena from the infrared channels of the Pressure 
Modulated Infrared Radiometer (PMIRR) and the Thermal Emis- 
sion Spectrometer (TES). These will supplementthe daytime global 
synoptic and high-resolution local visible images from the Mars 
Observer Camera (MOC). Similarly, the investigation of the global 
water vapor transport is also limited due to observational and orbital 
constraints. Mars Observer is more capable of yielding seasonal 
cycles of the martian atmosphere than the shorter-term atmospheric 
and surface phenomena that occur on a timescale of minutes and 
hours to days. These deficiencies arc addressed by anew Discovery- 
class mission that is currently being investigated. 

Mars Operational Environmental S atcllite (MOES) is a Discov- 
ery conccptmission that is designed to observe the global short-term 
weather phenomena on Mars in a systematic fashion. Even after the 
Mariner, Viking, and, soon, Mars Observer missions, crucial as- 
pects of the martian atmosphere will remain unobserved systemati- 
cally. Achieving a better understanding of the cycles of dust, water 
vapor, and ices on Mars requires detailed information about atmos- 
pheric transports of. those quantities associated with the weather 
systems, particularly those arising in mid latitudes during fall and 
winter. It also requires aquantitative understanding of the processes 
responsible for the onset and evolution of dust storms on all scales. 
Whereas on Earth the system of geosynchronous and polar orbiting 
satellites provides continuous coverage of the weather systems, on 
Mars the time history of important events such as regional and 
global dust storms remains unobserved. 


To begin with, we desire observations of Mars and its atmos- 
phere at a spatial resolution of about 3-5 km for visible and infrared 
imaging (for nighttime coverage) so that the evolution cycles of the 
dust storms and cloud systems on different scales can be observed 
both during the day and night In addition, we desire vertical profiles 
of temperature, water vapor, and dust globally, but at a lower 
horizontal resolution, somewhere between 30 and 75 km, corre- 
sponding to a grid point resolution between 1 ° and 2.5 °, comparable 
to the resolution of the numerical models. For the Earth atmospheric 
weather models the initial state is determined from the global 
network of surface and upper stations and augmented with the 
satellite observations of winds (cloud drift) and vertical temperature 
profiles at asynoptic times. On Mars, lacking a synoptic observation 
network of stations, all the data need to be determined from the 
satellite observations, which arc difficult to attain at the same 
observation time at all locations, which necessitates the use of 
assimilation of these observations in the numerical model. The 
min i m u m temporal resolution of these observations necessary to 
resolve the diurnal and the semidiurnal components adequately is 
about three hours (one eighth of the martian day, sol). If these arc to 
be obtained from a satellite in a low- altitude orbit as opposed to an 
aerosynchronous satellite, the temporal resolution in essence de- 
fines a maximum allowable altitude on the orbit of the single 
satellite. The altitude of the orbit also determines how much of the 
planet one is able to observe for different orbit inclinations. 

The MOES mission concept design is based upon these consid- 
erations and the payload is chosen to meet the observational require- 
ments. With its emphasis on short-term weather and time-of-day 
variations, MOES complements Mars Observer and Mars *94 mis- 
sions, which are optimal for seasonal-scale variations. A single 
MOES spacecraft deployed in a retrograde low-inclination (155°), 
216-minute orbit will yield observations of Mars surface and atmos- 
phere within ±75° latitudes up to 8 times per martian day. To 
observe the polar latitudes at the same time will require an addition- 
al spacecraft in a polar orbit that can be Sun-synchronous. Launch 
is on a Delta II 7925 and can take place in any of the Mars windows 
that recur. Orbit insertion at Mars is relatively easy due to the low 
inclination of the orbit, and no plane changes are required such as 
for MO. The candidate spacecraft is three- axis stabilized and is 
nadir oriented to provide the instrument coverage. 

The MOES primary scicncepayload consists of two instruments, 
Mars Imaging Radiometer (MIR), a multispectral imaging system 
(five channels), and Microwave Atmospheric Sounder (MAS), a 
multispectral microwave sounder to determine the vertical profiles 
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of temperature and water vapor by passive sounding. Both are 
identical twins of current Earth satellite counterparts with the 
exception of the filter set There are two’ instruments that currently 
provide the capability that MIR is required to have: the AVHRR 
being flown on the NO AA polar orbiting satellites and the VIRS on 
TRMM (to be launched in 1996). For Mars imaging applications we 
desire two visible (0.45 and 0.59 pm) and three infrared wave- 
lengths (6.9, 7.7, and 9.3 pm). The infrared channels are chosen for 
optimal coverage of the surface, water vapor, and atmospheric dust 

The Microwave Atmospheric Sounder is similar to the 1 83-GHz 
water vapor sounding instrument SSM/T2 currently in operation on 
the DMSP satellites and also has a separate receiver (230.5 GHz CO 
line) for deducing the martian atmospheric temperature profiles. 
Use of the CO channel for determining thermal structure is superior 
to the use of passive infrared techniques as it is unaffected by 
atmospheric dust and water vapor. The usefulness of the 1 83-GHz 
water vapor line to determine the atmospheric water vapor content 
on Mars has already been demonstrated by [1] from Earth-based 
observations using the VLA. A combined analysis of the infrared 
radiances and the temperature and moisture profiles will yield 
estimates of the dust content The primary payload is thus capable 
of satisfying the observational requirements. 

The visible and infrared imaging provides coverage of all mar- 
tian latitudes within ±75° up to eight times per sol by cross-track 
scanning with a nadir resolution of 4 x 4 km. Spatial coverage of 
temperature and moisture profiles is also obtained by cross-track 
scanning, but at a lower resolution of about 40 x 40 km at nadir. 
Additionally, MOES is occulted as seen from Earth on every orbit 
due to its low inclination, and provides independent estimates of 
high-vertical-resolution temperature profiles of the martian atmos- 
• phere that will be used for calibrating the microwave temperature 
profiles and dust abundance. 

Finally, the nominal orbit chosen for this mission will place 
MOES in sight of all the surface stations (within ±75°) planned to 
be deployed by MESUR and MARS NET for as much as 40 min per 
orbit, ormore than 5 hr per martian day, allowing ample opportunity 
to collect the data using its omnidirectional antennas for relay to the 
DSN using the High Gain Antenna. For coverage of the polar 
latitudes and also to receive the data from the polar surface stations, 
a second MOES spacecraft in a polar orbit is required and the two 
together would provide complete coverage of the planet Thus 
MOES mission(s) complements) not only the presently planned 
Mars orbiting missions in terms of science, but will also provide a 
global context for the MESUR observations, as well as be capable 
of providing a relay service for the surface network data. 

The mission operations plan is simple and essentially fixed: The 
two nadir-looking instruments observe continuously and the data 
are stored in onboard memory for transmission to the Deep Space 
Network (DSN) utilizing a single 34-m station pass per day. The 
HGA can support data rates from 587-17.3 kbps using the 
34-m DSN stations for Earth-Mars range between 0.45-2.62 AU. 
We envision routine spacecraft operations to be handled by the Jet 
Propulsion Laboratory (JPL) and the science data to be processed at 
the University of Wisconsin-Madison and made available in real 
time to the scientific community as resources. 

To su mm arize, MOES will yield spatially and temporally re- 
solved observations of atmospheric and surface temperatures, and 
vertical profiles of temperature and water vapor, as well as esti- 
mates of dust abundance at all longitudes. These observations will 


enable determination of the atmospheric transports by the martian 
general circulation, enable better characterization of the martian 
surface properties through infrared and visible observations, and 
provide crucial data for the initialization and verification of martian 
weather and climate models. This the scientific return is comple- 
mentary to the other missions to Mars of the past and those currently 
underway. By uti li z in g existing Earth applications instruments for 
martian application, substantial cost savings are expected to be 
realized. 
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been provided by TRW (C. Lillie) while mission design support has 
been provided by S. Weinstein (JPL). 
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A STUDY OF CLOUD MOTIONS ON MARS, H: SOME 
EXAMPLES OF OBSERVATIONS AND ANALYSES SINCE 
1969. L. J. Martin, Lowell Observatory, Flagstaff AZ 86001, 
USA 

The purpose of this paper is to remind everyone of the varied 
database on cloud motions that does exist but has often been under- 
used in dynamical studies. Some striking examples of these cloud 
motions will be described. As the title implies, an earlier study was 
made by the author, with W. A. Baum, to analyze observations of 
martian cloud motions [1]. That research was supported by the 
Mariner missions in the late 1960s under a contract with the Jet 
Propulsion Laboratory. Compared to data acquired since that study, 
the data used were weak and the conclusions were limited. This was 
primarily because the only observations available were sparse and 
noncontiguous. That project was completed in 1969; during that 
same year, Lowell Observatory, with NASA support, initiated a 
more comprehensive observing program under the direction of 
W. A Baum [2] . This International Planetary Patrol Program estab- 
lished a network of observing stations around the world that could 
provide nearly continuous observations of Mars throughout its 
apparitions. The primary objective of this network was to overcome 
the difficulties in attempting to document the motions of atmos- 
pheric activity using intermittent data. The network aspect of the 
Planetary Patrol tapered off after only three Mars oppositions, but 
observations from that period still provide the best data available, 
from any source, on martian cloud motions. 

Using the Planetary Patrol data from 1969, the hourly motions of 
the north polar hood clouds were mapped for each day for a period 
of several months [3]. During the same period, movements of the 
nonpolar clouds were also mapped hourly. Of greater interest, 
however, were the hourly cloud positions seen in the developing 
stages of the great planet-encircling dust storms of 1971 [4] and 
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1973 [5]. Maps of these motions were compiled for the first 20 days the half-hour that followed, this cloud dissipated, but motions of the 

of each of these storms. After 20 days, these storm clouds began to surrounding clouds had continued. Their apparent motions were 
merge into continuous dense hazes on these Earth-based images, radially outward from the location of the more discrete (but then 
This made it impractical to continue delineating individual clouds, gone) rillecloud. These images alsosuggestmovementof whatmay 

but because of those early studies, the motions of the 1971 and 1973 have been surface material, within the same 1-hrintervaL 

planet-encircling storms are better documented than the motions of A much smaller cloudmotion was seen on two consecudve high- 

any other martian dust storms. Daily diurnal expansions of these resolution Viking frames. The movement of the cloud and its shad- 

storms showed them moving farther west each day, although appar- ow within the 4.5-s interval between frames is probably another 

ently retreating to the east before rotating into view each morning, uniqueobservationforthe Vikingmission.Thiscloudalsomayhave 

During these early stages of storm development, condensate clouds been an eruption [10] or a dust devil, although it was so grp«11 that 

(blue-filter images) were often seen in association with, but usually it could not have been identified on any but the highest-resolution 

separate from, the dust clouds (on red-filter images). Viking images. 

The extent of transport of dust and condensates within these Other examples of cloud motion observed by Viking include the 

cloud motions is not straightforward. Since it is unlikely that dust diurnal “bore waves” over the Tharsis area [11,12], a local dust 
particles moved at rotation speeds, we assume that the dust clouds storm in the Chryse area [1 3], and the much larger regional dust 
brightened diurnaUy by a mechanism that was triggered by solar storm between 60° and 70°N that had been described as a moving 
heating or some type of thermal tides. Therefore, observed cloud “cold front” [14]. The Viking orbiters recorded a large number of 
motions are probably not a measure of motions of materials that different types of clouds, but only a small percentage of them on the 
make up the clouds. Likewise, vertical motions of cloud particles same day or consecutive days, giving us only sporadic glimpses of 
cannot be measured, although they are readily apparent on the more measurable motion. However, these data on moving clouds provide 
detailed Viking images. Because initial clouds of the 1973 encir- graphic examples of the nature of atmospheric transport on Mars, 
cling storm sprang up at different locations within a few hours or References: [1] Martin L. J. and Baum W. A. (1969) Final 

days of each other [5], and because there is evidence from several Report B, JPL Contract 951547, Lowell Observatory, Flagstaff, 
datasets that the martian atmosphere became dustier prior to the Arizona. [2] Baum W. A. et al. (1970) Icarus, 12, 437-439. 
first planet-encircling 1977 storm [6,7], it seems possible that the [3] Martin L. J. and McKinney W. M. (1974) Icarus, 23, 380-387. 
expansion of these storms feeds upon dust that is already in the [4] Martin L. J. (1974) Icarus, 22, 175-188. [5] Martin L. J. (1976) 
atmosphere rather than either horizontal transport or the continuous Icarus, 29, 363-380. [6] Martin L. J. et al. (1991 ) In LPI Tech. Rpt. 
uplifting of dust from the surface as the storm advances. Horizontal 92-02, 99-100. [7] Martin T. Z. and Richardson M. I. (1991 ) In LPI 
transport in itself probably requires too much speed and too much Tech. Rpt. 92-02, 101-102. [8] Martin L. J. and Zurek R. W. (1993) 
dust Continuous uplifting requires a path of continuous surface JGR, 98, 3221-3246. [9] Martin L. J. and Baum W. A. (1978) Bull 
locations that are both sources of dust and conducive to its uplift Am. Astron. Soc., 10, 551. [10] Martin L. J. (1981) Thirdlntema- 
Furthermore, given the special circumstances that are required to tional Colloquium on Mars, LPI Contribution 441, 145. 
initiate a storm (based upon their infrequency [8]), it seems unlikely [11] Hunt G. E. et al. (1980) Nature, 286, 352-364. 
that they could form initial clouds at widely separated locations (as [12] Kahn R. and Gierasch P. (1982) JGR, 87, 867-880. 
in 1973) withoutsomecommon conditions. Since surface conditions [13] James P.B. and Evans N. (1981) GRL, 8, 903-906. [14] James 
do not appear to be similar, it seems possible that the atmosphere’s P. B. (1983) In Recent Advances in Planetary Meteorology (G. E. 
increasing dust load might be. The similarity of the 1973 storm. Hunt, ed.), 85-99, Cambridge, 
observed from Earth, and the early 1977 storm, observed by Viking 

[9], allowed us to merge the data from sources with vastly different 

scales to provide a more comprehensive analysis of these events. EVIDENCE FOR DUST TRANSPORT IN VIKING IR 

The Viking data include several other examples of cloud motion, THERMAL MAPPER OPACITY DATA. T. Z. Martin, Jet 

including a few cases where imaging of the same cloud was repeated Propulsion Laboratory, CalifomialnstituteofTechnology, Pasadena 

during the same orbit providing hourly measures of motion. Most of CA91109, USA. 
these cloud events have been analyzed individually, but a compre- 
hensive catalog of all the cloud motions recorded by Viking has not Global maps of 9-pm dust opacity derived from radiometric 

yet been compiled. This is also true for recent Mars apparitions observations made by the Viking Qrbiter IR Thermal Mapper 

observed from Earth, although in the absence of a well-funded instruments have revealed a wealth of new information about the 

observing network, the number of documented cloudmotion events distribution of airborne dust over 1.36 Mars years from 1976-1979 

is limited. ^ . [1,2]. In particular, the ch anging dust distribution during major dust 

A unique Viking observation of a motion event was recorded storms is of interest since the data provide a point of contact with 
during a single orbit at half-hour intervals. Several clouds, identi- both Earth-based observations of storm growth and with global 
fiable by their shadows, were seen over the south flank of Arsia circulation models. 

Mons on the earliest of the three sequences of images. Half an hour The set of maps have a time resolution of 5° in L, (areocentric 

later, on the second sequence, a small, bright cloud appeared solar longitude) necessitated by the need to cover a large area in a 
directly over a rille and near the larger, earlier clouds (which showed singlemap; they are consequently not ideal for the detailed study of 
movement since the earlier sequence). Because of the sudden ap- dust transport. A different approach, re taining the full- time resolu- 
pearance andspecial location of this new cloud, I had suggested that tion of the date acquisition process, is to create a set of opacity time 
it might have been some type of eruption [10], although the possi- histories for selected sites of interest An example of one of these 
bility of something else like adust devil cannotberuled out Within histories, for a 10° x 10° area centered on latitude -40°, longitude 
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180°, is shown below for the L, range 200°-350°. Crosses indicate 
data from the instrument on Viking Orbiter 1; triangles are from 
V02. Although it is uncommon to obtain time resolution finer than 
about one day, and keeping in mind the effect of data coverage gaps 
that occur, several interesting points can be established with the set 
of opacity/time plots. There is a difference in the onset of the 1977b 
storm between northern and southern latitudes. Opacities remain 
high far longer at equatorial latitudes than in southern middle 
latitudes. Opacities at the Viking Lander sites were considerably 
smaller than those near the equator or to the south. 


- 40/180 
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Fig. 1. 

References: [1] Martin T. Z. (1986) Icarus , 66, 2. [2] Martin 
T. Z. (1993) JGR, in press. 


RETRIEVAL OF TEMPERATURE PROFILES FROM 
MARTIAN INFRARED SPECTRA. W. W. McMillan 1 ,! C. 
Pearl 2 , andB. J. Conrath 2 , 'NAS/NRCNAS A Goddard Space Flight 
Center, Code 693.2, Greenbelt MD 20771, USA, 2 NAS A Goddard 
Space Flight Center, Code 693.2, Greenbelt MD 20771, USA. 

Introduction: In theory, a thermal emission spectrum of Mars 
(200-1600 cm* 1 ) contains a wealth of information about the com- 
position and thermal structure of the atmosphere and surface. 
However, in practice, the ability to retrieve quantitative informa- 
tion from the spectrum is limited by a number of factors including 

(1) the ill-posed nature of the spectral inversion problem and the 
resulting nonuniqueness of all solutions; (2) assumptions built into 
any spectral inversion program; (3) uncertainties in surface pres- 
sure, temperature, and emissivity; and (4) uncertainties in the 
optical properties of atmospheric aerosols. Below, factors (1) and 

(2) are discussed as fundamental limitations on temperature re- 
trievals from Mariner 9 Infrared Interferometer Spectrometer (IRIS) 
spectra. 

In preparation for the Mars Observer mission and the return of 
tens of thousands of infrared spectra per day from the Thermal 
Emission Spectrometer (TES), we have developed a fast inversion 
algorithm to retrieve temperature structure and aerosol opacity from 


the infrared spectra. The derived atmospheric models will be used 
to provide the atmospheric contribution to the TES spectra, so that 
the thermal emission spectra of the underlying surface can be 
determined for making mineralogical identifications. As a test of 
our algorithm we are undertaking a systematic analysis of the entire 
Mariner 9 Infrared Interferometric Spectrometer (IRIS) dataset of 
21,000+ spectra. While portions of the IRIS dataset have been 
previously analyzed, the lack of a speedy and robust algorithm to 
invert the IRIS spectra to retrieve temperature profiles and aerosol 
opacities has been a major impediment to a truly comprehensive 
analysis. 

Technique: Generally, the most prominent feature of an IRIS 
spectrum is the 15-pm (667-cm* 1 )bandofC0 2 (seeFig. 1). Because 
of our understanding of the C0 2 molecule, its optical properties can 
be calculated and the 15-pm band inverted to obtain an estimate of 
the vertical thermal structure of the martian atmosphere. Indepen- 
dent temperature profiles can be obtained by separately inverting 
the two sides of the band (P and R branches). In the absence of any 
other atmospheric opacity, these profiles should be identical to 
within the rms error resulting from measurement noise propagation 
(1-2 K) and uncertainties in the C0 2 absorption properties. How- 
ever, in practice, additional wavenumber-dependent extinction due 
to absorption and scattering by atmospheric aerosols will cause the 
retrieved P and R branch profiles to differ. These differences will be 
exploited to assist in deriving information on opacity resulting from 
the presence of the particulates. 

W e use a constrained linear inversion (maximum entropy) tech- 
nique in an iterative scheme to solve the radiative transfer equation 
to retrieve P and R branch temperature profiles for each spectrum. 
Given a temperature profile, the transmittances required for the 
radiative transfer equation are interpolated in pressure, tempera- 
ture, and absorber mass from a lookup table of homogeneous path 
transmittances, using the Curtis-Godson approximation for the 
effective temperature and pressure. The homogeneous path trans- 
mittance lookup table is precomputed using an accurate line-by-line 
program. An iterative scheme is required because of the nonlinear 
relationship between temperature and measured radiances. For a 
spectrum taken through a clear atmosphere, a judicious choice of the 
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damping factor and a “good” initial guess lead to convergence in 
3-5 iterations. 

Temperature Retrievals: The retrieval of a temperature pro- 
file proceeds by first selecting a set of wavenumbers for use in the 
inversion. The radiances measured at these wavenumbers (IRIS 
spectrum) are used in the iterative scheme until the computed 
radiances for a retrieved profile are within the noise equivalent 
spectral radiance (NESR = 5 x 1(H W cm 2 sr cm* 1 ) of the IRIS 
spectrum. The wavenumbers are chosen to provide the maximum 
coverage in altitude with a uniform distribution of contribution 
functions. An appropriate set of wavenumbers and the respective 
contribution functions for the P-branch are shown in Fig. 2; a similar 
set of contribution functions is selected for the R branch. Retrieved 
P- and R-branch temperature profiles for the clear IRIS spectrum of 
Fig. 1 are shown in Fig. 3. 

At the 2.4-cnr 1 spectral resolution of the Mariner 9 IRIS, for 
nadir viewing, the center of the 15-pm band of C0 2 probes near the 
0.2-mb region of the martian atmosphere. This is evident from the 
broad peak of the contribution function for 668.45 cm- 1 (the IRIS 
sample nearest the band center) shown in Fig. 2. The contribution 



Fig. 2. Normalized contribution functions. 



Fig. 3. P- and R- branch retrieved temperature profiles. 
DAS = 13507967. 


function indicates the pressures sampled at a given wavenumberfor 
a given vertical temperature profile. Above the peak of the 668.45- 
cxn- 1 contribution function, the information content decreases as 
shown. Thus, the atmospheric temperature profile can only be 
accurately retrieved up to the vicinity of the pressure at the peak of 
the 668.4-cnr 1 contribution function; for a 10-km scale height this 
is roughly 35 km. 

Vertical structure that is significantly smaller than the half- 
width of the contribution functions (about a half scale height) is 
unresolvable. Thus, including additional wavenumbers to provide a 
denser vertical spacing of contribution functions will not substan- 
tially improve the resolution of a retrieved temperature profile. 

Our inversion technique is sensitive to the initial guess at the top 
of the atmosphere; we take the initial temperature profile above the 
peak of the 668.45-cm* 1 contribution function to be isothermal at the 
brightness temperature of the IRIS spectrum at 668.45-cm* 1 . We 
have found that such an initial guess minimizes the spectral residu- 
als at band center. The lower part of the initial temperature profile 
is estimated by assigning the brightness temperatures at selected 
wavenumbers to the pressures where the corresponding contribu- 
tion functions peak. 

Conclusion: Within the limitations discussed above, we have 
demonstrated the success of our inversion algorithm for the retrieval 
of temperature profiles for clear IRIS spectra like Fig. 1. Work 
continues on the modeling of aerosol extinction for the retrieval of 
temperatures and aerosol opacities from more dusty IRIS spectra. 
The simultaneous analysis of IRIS spectra of the same geographic 
location, but different viewing geometries, will be used in the 
characterization of the optical properties of airborne particulates. 

DUST TRANSPORT IN THE MARTIAN ATMOSPHERE. 
J. R. Murphy, NASA Ames Research Center, Moffett Field CA 
94035-1000, USA. 

Dust in suspension within the martian atmosphere is an impor- 
tant driver of the atmospheric thermal and, ultimately, dynamical 
states. By virtue of its presence in the atmosphere, this suspended 
dust is susceptible to transport by the winds and thus the location 
from which dust is lifted from the surface is not necessarily the 
location at which it will return to the surface. Such surface dust 
lifting and subsequent transport, when accumulated over daily, 
seasonal and annual time intervals define the martian dust cycle [1], 
The atmospheric transport element (and equally surface lifting and 
redeposition) of this cycle is strongly dependent upon season and the 
quantity of dust in suspension. The seasonal dependence arises due 
to the dominance of particular components of the atmospheric 
circulation atparticular locations (latitudes) at various times through 
the year. The dependence on dust abundance is due to the amount 
of dust available for transport, but it is also due to the influence that 
the dust has on the intensity of the circulation. As dust abundance 
and therefore its radiative influence increases, some components of 
the circulation intensify (Hadley circulation, thermal tides), while 
some components may decline in intensity (baroclinic waves, con- 
densation flow) [2,3]. For these reasons, our ability to understand 
the martian dust cycle is dependent on our ability to define the 
contributions that various circulation components, under varying 
seasonal and dust-radiative forcings, play in the transport of dust 
around the planet 
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Below we will discuss a number of components of the atmos- 
pheric circulation and what their impacts might be upon net dust 
transport The points made here are based in part on observations, 
on numerical modeling, and on extrapolations of our understanding 
of transport processes in the terrestrial atmosphere. The presenta- 
tion here is less quantitative than one might prefer, but it helps to 
frame the questions that need to be addressed to produce the 
quantification we desire. 

The best-characterized transport mechanism, and possibly most 
important in a global sense, is provided by the Hadley circulation. 
This thermally direct cell, with its ascending branch situated near 
tiie subsolar latitude at any particular time of year, can span a wide 
latitudinal range and provide an efficient mechanism by which dust 
lifted from a confined area can be redistributed over a large fraction 
of the globe. This is seemingly what occurs during global dust 
storms. During these events, which preferably begin in the southern 
hemisphere, during southern spring and summer, dust fills the 
atmosphere to a depth of several scale heights and covers most of the 
globe. The cross-equatorial Hadley cell responds by intensifying 
and expanding meridionally, with dust raised at southern subtropi- 
cal and middle latitudes spreading into the northern hemisphere at 
elevations one to several scale heights above the surface. This 
intensified circulation, which rapidly transports dust, is also prob- 
ably involved in the raising of dust from the surface. Numerical 
modeling has shown the Hadley circulation to be quite efficient at 
transporting dust from a prescribed southern hemisphere subtropi- 
cal/middle latitude source to northern middle latitudes in the span 
of only several sols [4,5]. This rate of dust spreading is consistent 
with observations. 

Eastward-migrating mid- and high-latitude waves are present 
during autumn, winter, and spring. These baroclinic waves arise in 
response to the thermal gradient between the warm equator and cold 
pole at these seasons. These waves might play an important role in 
the transport of dust to the polar regions. The airflow within these 
waves exhibits poleward flow of warm air in advance of the east- 
ward migrating wave, and the equatorward flow of cold air in the 
wave’s lee. Since the air equatorward of the seasonal cap should in 
general possess a greater dust content than air over the cap, since the 
C0 2 cover prevents the lifting of dust, one might expect baroclinic 
waves to provide a net poleward transport of dust. Since the 
baroclinic waves have a finite latitudinal extent, as the seasonal cap 
grows, the northward extent to which the baroclinic waves can 
transport dust will migrate equatorward. Similarly, as the seasonal 
cap retreats in the spring, the poleward limit of the dust transport 
will advance toward the pole. In this way, baroclinic waves can play 
a major role in the incorporation of dust into the seasonal cap. They 
also can play a role in transporting dust to the pole for possible 
inclusion in the polar laminae. This dust can either come directly 
from much lower latitudes, or it can be dust that was previously 
incorporated into the seasonal cap earlier in the season, which then 
becomes susecptible to resuspension and transport after the sea- 
sonal cap has sublimed away. 

Viking Lander 2 meteorological and visible opacity observations 
at a northern mid-latitude site (47°N) do suggest a correlation 
between the passage of baroclinic waves and variations in the 
overhead opacity [6] . Opacity increases during the poleward flow in 
advance of the wave trough, and after the trough passage opacities 
decrease in the equatorward flow behind the trough [7], as our 
conceptual model suggests. Numerical modeling of dust transport in 


the martian atmosphere also shows such correlations (Fig. 1). In 
fact, baroclinic waves are the most important element for transport- 
ing dust to high northern latitudes during a simulated global dust 
storm with a specified southern subtropical dust source [8]. 

Baroclinic waves exist in years with and without global dust 
storms, and thus their effects upon net dust transport to high lati- 
tudes might well be more important than transport that might arise 
from the episodic global dust storms (discussed further below). 
Numerical modeling of dynamical processes has pointed out that 
baroclinic wave activity is more intense in the northern hemisphere 
than in the southern hemisphere [2]. These differing amplitudes 
might be a response to the model topography (which may be incor- 
rect!) differences in the two hemispheres. However, another reason 
for hemispheric differences is probably due to the greater latitudinal 
extent of the southern hemisphere seasonal C0 2 cap and concurrent 
changes in the zonally averaged thermal and wind fields, which 
reduce the baroclinic instability [2]. These hemispheric differences 
in wave amplitude should then result in hemispheric differences in 
wave transport of dust to high latitudes. The persistence of a south 
polar residual C0 2 cap through the summer months during some 
years could possibly be a manifestation of this difference, dust 
contamination of the north cap producing a lower albedo, that 
prevents retention of a C0 2 cap there. 

A component of the martian atmospheric circulation that distin- 
guishes it from its terrestrial counterpartis the substantial seasonal 
net mass transfer between the polar regions. This “condensation” 
flow arises due to the condensation upon the winter pole of a 
significant fraction (-30%) of the atmospheric mass [9,10]. This 
flow is manifested by a net mass -averaged meridional flow from the 
hemisphere of the subliming seasonal cap to the hemisphere of the 
condensing seasonal cap. The hemispheric difference in the sizes of 
the seasonal caps, caused by Mars* large orbital eccentricity, im- 
plies differences in the magnitude of the condensation flow with 
season. Such differences, in conjunction with seasonal variations of 
dust loading during years with and also without global-scale storms, 
might lead to a net transfer of dust from one hemisphere to another. 
The condensation flow magnitude during southern autumn/winteris 
greater than the magnitude during northern autumn/winter. The 
atmospheric dust loading is generally greater during northern au- 
tumn/winter than during the same southern seasons. This leaves 
open the question of net transport, since the correlation between 
flow magnitude and dust loading variations might balance out over 
an annual cycle. 



Fig. 1. The model-calculated near-surface meridional wind and column- 
integrated visible dust opacity at a grid point corresponding to the Viking 
Lander 2 northern mid-latitude site. The positive cocrelationbetween the variables 
represents a net poleward transport of dust by low-frequency baroclinic waves. 
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Dust transport induced by atmospheric tidal motions is very In this presentation I have attempted, albeit rather crudely, to 
poorly defined. There are several scenarios in which dust transport describe various components of the martian atmospheric circulation 
by the tides can be envisioned. During the developing stage of a and their implications for dust transport The seasonality of various 
global dust storm, tidal amplitudes rapidly intensify. This amplifi- components must surely be considered when attempting to define 
cation occurs during a time period in which the atmospheric dust net annual and multi-annual effects of individual components, as 
load is rapidly increasing due to sources located at southern subtrop- should the seasonality of suspended dust abundance. Quantification 
ical latitudes. The large meridional (and zonal) flows that accom- of these various processes, both from numerical models and the 
pany the tidal amplifications can produce large dust fluxes if dust is Mars Observer dataset will tremendously expand our understand- 
presentatthc appropriate vertical levels. Time-constant tidal ampli- ing of the martian dust cycle and the relative importance of its 
tudes are less likely to produce a net transport of dust; however, dust various constituent parts. 

particle sedimentation, in conjunction with nonamplifying tidal References: [1] Kahn et al. (1993) Mars, 29, Univ. Arizona, 
flow, can produce net dust transport It has been shown [11] that [2] Haberle ct al. (1993) JGR, 98, 3093-3123. [3] Barnes ct al. 

inclusion of tidal forcing in a zonal-mean model during a simulated (1993) JGR, 98, 3125-3148. [4] Haberle et al. (1982) Icarus, 50, 

global duststorm(largedustloading)enablestheHadley circulation 322-367. [5] Murphy ctal. (1993)76/?, 98, 3197-3220. [6] Tilman 

to extend to higher latitudes in the winter hemisphere, with a etal. (1979) JGR, 82, 2947-2955. [7] Ryan and Henry (1979) JGR, 

concurrent shift of the northward extent of dust transport However, 84, 2821-2829. [8] Murphy etal. (1993) in preparation. [9] Pollack 

this expansion does not extend into polar regions. The rapid inten- et al. (1990) JGR, 95, 1447-1474. [10] Pollack et al. (1993) JGR, 
sification of tidal amplitudes at the onset of dust storms and the very 98, 3149-3181. [llJZurek and Haberle (1988) JAS, 48, 1005-1023. 

large wind speeds thatresult are certainly an inviting mechanism to [12] Zurek and Christensen (1990 )Bull AAS, 22, 1075. [13] Barnes 

invoke to explain the rapid growth of the observed storms. Much (1990) JGR, 95, 1381-1400. [14] Jakosky andMartin (1987) Icarus, 
more work is required to better define the importance of such a 72, 528-534. 
process. 

Regional types of flows can also produce net dust transport Eas t- 

west circulations at northern subtropical latitudes have been in- THREE-DIMENSIONAL NUMERICAL SIMULATION OF 
voked to explain an inferred removal of dust from high-thermal- NEAR-SURFACE FLOWS OVER THE MARTIAN NORTH 
inertia, low-elevation, dark areas, and subsequent redeposition of POLAR CAP. T. R. Parish 1 and A. D. Howard 2 , 1 Departmentof 
this same dust upon low-thermal-inertia, higher-elevation, visibly Atmospheric Science, University of Wyoming, Laramie WY 82070, 
brighter locations [12]. Such cast-west, or Walker type, circulations USA, department of Environmental Sciences, University of 
could themselves be driven by the heating contrasts engendered by Virginia, Charlottesville VA 22903, USA. 
the surface properties outlined above. 

Diumally forced upslope/downslope flows along the numerous Introduction: Measurements made by Viking Lander VL-2 
sloping regions on Mars are another local orregional factor that can (48°N) have shown that the near-surface wind and temperature 

be invoked to produce a net dust transport The Hellas Basin, regime on Mars displays striking similarities to terrestrial counter- 
characterized by sloping faces on the order of 1000 km in length, parts [1]. The diurnal radiative cycle is responsible for establish- 
have been observed as preferential locations for the appearance of ment of a well-defined thermal circulation in which downslope 
suspended dust The southeast flank of Tharsis is another location (katabatic) flows prevail during the nighttime hours and weak 
of apparently preferred dust raising. The factors involved in and the upslope (anabatic) conditions prevail during the daytime. Previous 

importance of such relatively small-scale (mesoscale) processes work [2] has indicated that the slope flows are much like those found 

upon the martian dust cycle have not been thcoreticaly dealt with at on Earth, particularly the katabatic winds, which show striking 

similarities to drainage flows observed over Antarctica. The low- 

A final class of motions that can be invoked as transport-capable level wind regime appears to be an important factor in the scouring 
is what I will term “special circumstances.” These features are those of the martian landscape [3]. The north polar cap shows evidence of 
that might be sporadic in nature but can possibly have a very large eolian features such as dunes, frost streaks, and grooves from Viking 
impact upon net transports over a range of timescale. The best imagery. The direction of the prevailing wind can in cases be 
example of such a process is the circumstances surrounding a win ter inferred from the eolian features. 

polar warming similar to that observed during the second 1977 Here we examine the thermally induced flows that result from 
global dust storm. As [13] has pointed out, whatever dynamical the radiative heating and cooling of the martian north polar region 
process is fundamentally responsible for the warming, it must using a comprehensive three-dimensional atmospheric mesoscale 
involve poleward flow at high altitude and descent over the pole numerical model. The same model has been used previously for 
(which adiabatically generates the observed wanning). The persis- simulation of Antarctic katabatic winds [4], The model equations 
tence of the warming, -30 sols [14], and its intensity suggest that are written in terrain-following a coordinates to allow for irregular 
substantial meridional transport of various species (dust, water) terrain [5]; prognostic equations include the flux forms of the 
from low latitudes to the polar region could occur during such an horizontalmomentum equations, temperature, continuity. A surface 
event This dust could then become incorporated into the polar cap energy budget equation is also incorporated in which the surface 
due to either direct sedimentation upon the surface or its incorpora- temperature is determined. Explicit parameterization of both teires- 
don into water or C0 2 “snow flakes," which subsequently fall to the trial (longwave) and solar (shortwave) radiarion is included. Turbu- 
surface. The viability of this transport process depends strongly lent transfer of heat and momentum in the martian atmosphere is 
upon the vertical and latitudinal distribution of dust at low latitudes assumed to follow the similarity expressions in the surface boundary 
during the time period of the warming. layer on Earth [6-8]. The terrain heights for the martian north polar 
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regionhavebeenobtainedfromthcU.S.GeologicalSurveymapand Effect of solar heating of the martian north polar region, 
digitized onto a 57 x 57 grid with a spacing of 75 km. The resulting Numerical simulations have also been conducted incorporating 

terrain map is shown in Fig. 1 . The vertical grid consists of 15 levels solar insolation over the irregular terrain over the northern martian 

(o = 0.998, 0.99, 0.98, 0.97, 0.96, 0.94, 0.92, 0.90, 0.85, 0.775, latitudes. The model terrain heights are as shown in Fig. 1. Model 

0.70, 0.60, 0.50, 0.30, 0.10). The high resolution in the lower equations have been integrated for a period of three martian days to 

atmosphere is necessary to capture details of the boundary layer ensure that the wind and temperature regimes have had ample time 

flows. The lowest level corresponds to a height of approximately to adjust to the cycle of solar heating (solar constant is set to 600 W 

20 m above the ground, the second level 100 m. m~ 2 ). It was found that the second- and third-day results of the model 

Model Results: Katabatic wind simulations over the north integration were nearly identical. Here the third-day results of one 

polar cap . To examine the intensity and pattern of drainage flows 
over the martian north polar region, a series of numerical experi- 
ments have been conducted. To isolate the drainage flows, the 
model is initialized about a state of rest. No ambient horizontal 
pressure gradients are assumed to exist in the model atmosphere at 
the start The model simulation is assumed to be representative of 
polar night conditions; solar radiation is neglected throughout the 
model domain. The motions that develop are therefore the result of 
the longwave radiative flux divergence at the surface and in the low 
levels of the atmosphere. The model equations are integrated for a 
24-martian-hr period, by which time near-steady drainage flows 
have developed. 

Figure 2 depicts the streamlines and wind speeds of the katabatic 
drainage at the end of the 24-hr integration period at the second a 
level (approximately 100 m above the surface). Note that the drain- 
age flows move radially outward from the various ridges atop the 
north polar region, being deflected some 40°- 60° to the right of the 
fell line due to the Coriolis influences. Confluence zones, regions 
where the drainage streamlines converge, can be identified as well. 

In such areas cold, negatively buoyant air collects; katabatic wind 
can be significantly enhanced along the confluence axis. The wind 
speeds of the drainage flows are comparable to what would be 
observed over Antarctica. The strongest wind speeds appear to be 
related to the steepness of the terrain. Confluence effects accountfor Fig. 1. Terrain contours (200-m increment) near the martian north polar 
only a minor enhancement in the simulated katabatic intensities. region adapted from U.S.G.S. map. 
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Fig. 2. Model simulations of streamlines (left) and wind speeds (right) of katabatic flows at the second a level after 24-hr integration period. 
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such experiment will be presented. In this experiment, the solar 
declination is assumed to be 25°, representative of midsummer 
conditions. It is assumed that the atmosphere is dust free. As before, 
the model is initialized about a state of rest to isolate the thermally 
induced slope flows in the lower atmosphere. Values of albedo and 
thermal inertia were taken from Paige [9]. 

Figure 3 illustrates the solar insolation received at the surface 
and corresponding surface temperatures. Local noon at this time is 
along the 60° meridian in Fig. 1. Note that significant spatial 
variations in the intensity of solar radiation reaching the surface are 
foundnear thepolarcapin response to albedodifferenccs.Tempera- 
tures range from 260 K along the southern edge of the 60° meridian 


in the model domain to around 190 Kin the nighttime section of the 
model domain situated aboutthe210° meridian. Examination of the 
surface temperature regime shows little diurnal variation over the 
north polar ice cap; the largest diurnal changes are seen in the 
midlatitude portion of the model domain, as expected. The temper- 
atures at the first a level (not shown) vary from approximately 
238 K along the 30° meridian to approximately 198 K in the 
nighttime sector of the model domain. 

The thermally induced slope flows become significantly modi- 
fied from the pure katabatic drainage pattern shown in Fig. 2. 
Anabatic motions become established in response to the intense 
solar heating of the sloping terrain. Figure 4 illustrates the stream- 



Flg. 4. Model simulations of streamlines (left) and wind speeds (right) of slope flows at the second O’ level for diurnal cycle experiment. 
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lines and wind speeds at the second o level for the time of the local 
solar forcing shown in Fig. 3. Note that especially along the 30° 
meridian, the streamline patterns are dramatically altered by the 
solar insolation as compared to Fig. 2. Upslope wind conditions are 
simulated at middle to high latitudes from approximately 3 1 5°-90°. 
As seen in the previous simulation, the strongest winds are associ- 
ated with the steepest terrain. Upslope flows reach 7 ms- 1 along the 
60° meridian. The intensity of the upslope flows appears to be 
weaker than their nocturnal katabatic counterparts. This is in agree- 
ment with numerous terrestrial observations such as over the Ant- 
arctic ice sheet. Note that the katabatic wind speeds simulated over 
the nighttime sector centered about 210° are only slightly weaker 
than the katabatic winds seen in Fig. 2, emphasizing the rapid 
response of the wind field to the radiative-induced temperature 
changes near the surface. 

Summary: Numerical simulations of the martian near-sur- 
face wind regime using a three-dimensional atmospheric model 
have shown that katabatic winds prevail over the sloping terrain 
during nighttime hours. Such winds are forced by the strong radia- 
tive cooling of the sloping surface, which sets up a downslope- 
directed horizontal pressure gradient force. Both the magnitude and 
direction of these flows are similar to terrestrial counterparts such 
as those over Antarctica. Introduction of solar radiation strongly 
perturbs the slope flows; anabatic conditions develop in middle to 
high latitudes during the daytime hours due to the solar heating of 
the sloping terrain. Wind speeds for the upslope flows are somewhat 
weaker than for corresponding katabatic flows. There appears to be 
a rapid transition from the katabatic to the anabatic flow regimes, 
emphasizing the primary importance of radiative exchanges at the 
surface in specifying the horizontal pressure gradient force. 

References: [ 1 ] Hess S. L. et al. ( 1 977) JGR , 82, 4559-4574. 
[2] Parish T. R, and Howard A. D. ( 1 993 )JGR, submitted. [3] Howard 
A. D. (1981) NASA TM 82385, 333-335. [4] Parish T.R. and Waight 
K. T. ( mi)Moru Wea. Rev., 1 15, 22 1 4-2226. [5] Anthes R. A. and 
Warner T. T. (1978) Mon. Wea . Rev., 106, 1045-1078. [6] Brost 
R. A. and Wyngaard J. C. (1978)7. Armor. Sci, 35, 1427-1440. [7] 
Busch N. E. et al. ( 1 976)7. Appl. Meteor., 15, 909-9 1 9. [8] Businger 
J. A. et al. ( 197 1)7. Atmos. Sci., 28, 181-189. [9] Paige D. A. (1992) 
LPSXX1II, 1013. 


THE ROLE OF ATMOSPHERIC HEAT TRANSPORT IN 
THE SEASONAL CARBON DIOXIDE CYCLE. J. B. 
Pollack 1 , R. M. Haberle 1 , J. R. Murphy 1 , J. Schaeffer 2 , and H. Lee 2 , 
! NASA Ames Research Center, Moffett Field CA 94035-1000, 
USA, Sterling Software, Inc., Palo Alto CA 94303, USA. 

We have carried out numerical experiments with a general 
circulation model (GCM) and energy balance model of the martian 
atmosphere to define the importance of heat transported to the polar 
regions in determining the amount of C0 2 condensed on the surface 
during the fall and winter seasons and the amount sublimated during 
the spring and summer seasons [1,2]. In so doing, we performed both 
sensitivity experiments, in which the dust opacity was varied over 
the full range of its observed values [3], and annual simulations, in 
which the dust opacity varied continuously with seasonal data, in 
accord with measurements taken at the Viking landers [3], 

Dust opacity represents the key variable for determining the 
contribution of atmospheric heat advection to the energy budget in 


the polar regions. The amount of heat advected to the winter polar 
regions increases monotonically as the dust opacity at low and 
middle latitudes increases. However, the increase is sharpest be- 
tween optical depths of 0 and 1 and tends to level off at still higher 
optical depths. Heat advection is more important at times of C0 2 
condensation than C0 2 sublimation, since the temperature gradi- 
ents are much steeper in the winter hemisphere than in the summer 
hemisphere. Because dust opacity is much higher during northern 
winter than during southern winter, atmospheric heat advection 
reduces the amount of C0 2 that condenses in the north by a much 
larger factor than it does in the south. 

Three wind components play key roles in transporting heat to the 
polar regions: the zonally mean overturning circulation (Hadley 
cell), the transient eddies (baroclinic eddies), and the condensation 
flow. As the dust opacity increases, the Hadley circulation strongly 
expands poleward (but its poleward expansion is ultimately limited 
by angular momentum considerations); eddy transport shifts 
poleward; and the condensation flow weakens. 

The importance of atmospheric heat advection can be illustrated 
by comparing an annual simulation that includes it and whose 
parameters have been adjusted to match the seasonal pressure 
variations measured at the Viking landers [4] with an identical 
simulation, in which atmospheric heat advection is ignored. Pres- 
sure differences of as much as 0.5 mbar are found between the two 
runs, i.e., heat advection is responsible for several tens of percent of 
the observed signal. The seasonal pressure curves show a greater 
degree of symmetry with heat advection removed and an even 
greater degree of symmetry when seasonal hemispherical shifts in 
the global atmospheric pressure distribution is ignored. 

References: [1] Pollack J. B. et al. (1990) JGR, 95, 1447- 
1 473. [2] Pollack J. B. et al. ( 1 993)JGR, 98, 3 1 49-3 1 8 1 . [3] Colburn 
D.S. et al. (1989) Icarus, 79, 159-189. [4] Tillman J. E. (1988) 
JGR, 93, 9433-9451. 


EDDY MIXING COEFFICIENT UPPER LIMIT DE- 
RIVED FROM THE PHOTOCHEMICAL BALANCE 
OF 0 2 . J. Rosenqvist 1 and E. Chassefiere 2 , 1 DESPA, Observatoire 
de Paris-Meudon, 92195 Meudon, France, 2 C.N.R.S., Service 
d* Aeronomie, B.P. No. 3,91371 Verrieres-le-Buisson Cedex, France. 

The Model: This work is based on the study of the photo- 
chemical balance of molecular oxygen in the martian atmosphere by 
using a one-dimensional model of photochemical reactions involv- 
ing species derived from C0 2 and H 2 0. 

J, C0 2 + h v CO + O 
J 2 0 2 + h v — > 2 O 

J3 H 2 0 + h v — » OH + H 

J 4 H 2 0 2 + h v 2 OH 
J 5 0 3 + h u O 2 0A) + O OD) 

J 6 H0 2 + h v — » OH + O 

R, C0 2 + 0( 1 D)^C0 2 + 0 
R 2 OOD) + H 2 0 -» 2 OH 
R 3 Q('D) + U 2 ->0H + H 
R 4 H + 0 2 + C0 2 -» H0 2 + C0 2 
R5 0 + 0 + C0 2 — > 0 2 + C0 2 

Rg O + 0 2 + C0 2 — > O3 + C0 2 

r 7 o+ho 2 ->oh + o 2 
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Rg 0 3 + H — ) OH + O 2 
Rg 0 + 0H->H + 0 2 
R 10 CO + OH -> C0 2 + H 
R n H + H0 2 — > 2 OH 
R 12 H + ho 2 — > H 2 + o 2 
R 13 H + H0 2 H 2 0 + O 

R]4 2 H0 2 — > H 2 0 2 + 0 2 

r 15 ho 2 + oh->h 2 o + o 2 
r 16 h 2 o 2 + oh -> ho 2 + h 2 o 

The mcxiel is basically similar to the one used by [ 1 ] for the study 
of the regulation of CO on a global scale, but the chemical rates are 
taken from [2] . In the present scheme, the regulation of molecular 
oxygen is studied over timescales of the order of its photochemical 
lifetime («30 yr), which is much shorter than typical escape timescales 
«0 2 ]/<t>(0) * 10 5 yr). Thus, the escape fluxes are fixed to the values 
given by [3] and [4]. We examine the calculated equilibrium abun- 
dances of 0 2 for given thermal, eddy diffusion coefficients and H 2 0 
profiles. The thermal profile is taken from [5] in the lower atmo- 
sphere. At higher levels, in order to include the diurnal and seasonal 
thermal profile variability, we have also used the IRTM data [6]. 
The pressure at the surface has been taken from [7] . In order to study 
the influence of both temperature and pressure profiles on the 0 2 
mixing ratio, we have made several tests corresponding to different 
martian seasons. The results show that the influence of pressure 
and temperature is quantitatively weak compared to the one of 
K and of the water vapor density [H 2 0]. Thus, in the follow- 
ing we have fixed the pressure at the surface to a value of 
7 mbar and we have used a unique standard thermal profile corre- 
sponding to a profile roughly averaged over the year, the season, and 
the day: T = 205 K at 0 km altitude, 175 K at 25 km, and 145 K at 
50 km. 

Variation of [O 2 ] with K: Increasing K will accelerate the 

transport of O atoms, large amounts of which are formed above 
40 km by the C0 2 photolysis, toward the lower atmosphere where, 
in turn, they will be destroyed by reacting with odd hydrogen 
radicals H0 2 (reaction R 7 ) and OH (reaction R 9 ). By taking into 
account only the dominant terms, the residual production rate of 0 2 , 
when subtracting balance equations of 0 3 and H0 2 , may be written 

P Q2 = k 5 [0] 2 [C0 2 ] + k,[0][0H] (1) 

Reactions R 5 and R 9 are the only terms entering the production 
of 0 2 averaged over a time larger than the 0 3 and H0 2 lifetimes. 
They work mainly above 40 km, R 5 being maximum between 60 and 
70 km, and R 9 at 50 km. As K increases, R 9 becomes more important 
relatively to R 5 . 

At altitudes greater than 60 km, the balance of O shows that 
oxygen atoms reform 0 2 . But when K is large, the O atoms are less 
abundant at these levels because some of them are transported down 
to lower altitudes. For a fixed water vapor vertical distribution, 
when K is large, the downward transport of O atoms is important. 
In this case, the production of molecular oxygen is inefficient. 
Conversely, production of 0 2 by the action of O with itself becomes 
important for small values of K (^lO 6 cm 2 s -1 )* However, for 
extremely low values of K (« 1 0 5 cm 2 s _1 ), the production of 0 2 again 
decreases. This fact can be explained by the very slow transport of 
O atoms produced by photodissociation in the upper atmosphere 
towards the 50-70-km altitude range where R 5 and R 9 occur. Thus, 


for a given distribution of H 2 0, the 0 2 mixing ratio is globally 
anticorrelated with the eddy diffusion coefficient when the latter is 
somewhat larger than 10 6 cm 2 s _1 . 

Variation of [O 2 ] with [H 2 0]: If the water vapor amount is 

large, the formation of odd hydrogen (H, H0 2 , OH, H 2 0 2 ) radicals 
is also large. Then, in order to ensure the balance of odd hydrogen 
radicals, the destruction of H* should increase. It shows that H 
atoms and H0 2 densities must globally increase. At this step, the 
equilibrium of H atoms over the whole atmosphere can be written 
as 

J k 10 [CO][OH]dz + J k,[o][OH]dz = 

J k 4 [0 2 ][H][C0 2 ]dz + J k 8 [0 3 ][H]dz 

Quantitative comparison of the first and second terms of the left side 
of the above equation shows that 

J k 10 [CO][OH]dz » J k,[0][0H]dz (3) 

whatever the input values of K and [H 2 0]/[C0 2 ]. In the same 
manner, the comparison of the right side of this equation leads to 

J k 4 [0 2 ][H][C0 2 ]dz » J k 8 [0 3 ][H]dz (4) 

Moreover, the first part of the left side of the equation is constant as 
shown by [1]. Consequently, the product [II] [0 2 ] is expected to be 
roughly constant. Thus, as the amount of water vapor increases, [0 2 ] 
decreases. 

Results: As has been mentioned previously by several authors 

[8-1 1], an interesting property of 0 2 is that its abundance is rather 
sensitive to the eddy diffusion coefficient K and the amount of water 
vapor in the atmosphere. Thus, the knowledge of the 0 2 abundance 
(1.1 -1.3 ICE 3 ) measured from Earth [12-14], allows an interesting 
constraint on K and H 2 0. This model being a stationary one- 
dimensional model, these constraints correspond to eddy diffusion 
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Fig. 1. Calculated 0 2 mixing ratio as a function of K and the total water vapor 
amount in the atmosphere. The scale indicates the ranges of 0 2 mixing ratios. 
The clear gray shaded area corresponds to expected values of the 0 2 mixing ratio 
(1.1 x ICE 3 - 1.3 xl(E 3 ). The white area corresponds to an 0 2 mixing ratio larger 
than 1.5 x 10~ 3 and the dark one to a 0 2 mixing ratio lower than 9 x 1(H. 
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sional model, these constraints correspond to eddy diffusion coeffi- 
cients and water vapor amounts averaged over time (several de- 
cades) and over space (whole planet). Figure 1 shows the range of 
the molecular oxygen mixing ratio as a function of the eddy diffusion 
coefficient and water vapor total amount The clear gray area corre- 
sponds to the 0 2 mixing ratios, which are consistent with the 
observations (0 2 = 1.1-1 .3 ' ICh 3 ). An upper limit of 9/pm for H 2 0 
may be derived from this figure. This limit (which is an average over 
the night and day) is in agreement with expected abundances of H 2 0 
[15]. On the other hand, from Fig. 1, the absolute upper limit on the 
eddy mixing coefficient is 3 ' 10 7 cm 2 s* 1 in the lower atmosphere 
(20-50 km). This value is consitent with recent estimates of K 
derived from observations [16,17]. If the water vapor amount is 
lower than 4 pr pm, K should be larger than 10 7 cm 2 s Conversely, 
if the water vapor amount is larger than 6/pm, K should be lower 
than 10 7 cm 2 s* 1 . This last result shows that, if compared with the 
observations [16,17], the mean diurnal H 2 0 amount must be larger 
than 4/pm. 
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DIAGNOSTIC CALCULATIONS OF THE CIRCULATION 
IN THE MARTIAN ATMOSPHERE. M. Santee and D. Crisp, 
Earth and Space Sciences Division, Jet Propulsion Laboratory, 
California Institute of Technology, Pasadena CA 91109, USA. 

The circulation of the martian atmosphere during late southern 
summer is inferred from observed atmospheric temperature and 
dust distributions. We use global maps of temperature and dust 
optical depth (-0-60 km) retrieved from a subset of the Mariner 9 
IRIS thermal emission spectra spanning L, = 343°-348° [1]. This 
thermal structure is characterized by a reversed meridional tem- 
perature gradient at altitudes above about 40 km, and temperatures 
that decrease from equator to pole at lower altitudes. Zonal-mean 
zonal winds are derived from the zonally averaged temperatures 
assuming gradient wind balance and zero surface zonal wind (Fig . 1 ). 
Both hemispheres have intense midlatitude westerly jets with ve- 
locities of 80-90 m s -1 near 50 km; in the southern tropics the winds 
are easterly with velocities of 40 m s- 1 near 50 km. The north-south 
atmospheric transport includes contributions from both the zonal- 
mean meridional circulation and large-scale waves. 

Their net effect can be approximated by the diabadc circulation, 
which can be derived from the atmospheric thermal structure and 


the net radiative heating rates. A radiative transfer model [2] that 
accounts for absorption, emission, and multiple scattering by par- 
ticles and nongray gases is used to compute the solar heating and 
thermal cooling rates from diurnal averages of the retrieved IRIS 
temperature and dust distributions (Fig. 2). At pressures below 
4 mbar, there are large net radiative heating rates (up to 8 K <H) in 
the equatorial region and large net radiative cooling rates (up to 
20 K d- 1 ) in the polar regions. These net radiative heating rates are 
used in a diagnostic stream function model that solves for the 
meridional and vertical components of the diabatic circulation si- 
multaneously. We find a two-cell circulation, with rising motion 
over the equator, poleward flow in both hemispheres, s inkin g mo- 
tion over both polar regions, and return flow in the lowest atmo- 
spheric levels (Fig. 3). The maximum meridional velocity is 4 m sr 
1 in the tropics at -55 km altitude, and the ma ximum vertical 
velocity is 4 cm s- 1 downward over the north pole at about 60 km 
altitude. This circulation is sufficiently vigorous that the meridional 
transport timescale is about 28 days. The assumption that the 
advective and diffusive transport timescales should be roughly 
comparable leads to a vertical eddy diffusion coefficient of about 30 



Fig. 1. Zonal-mean zonal winds, in ms -1 , calculatedfromthegradientthermal 
wind equation. Positive values represent eastward winds, negative contours are 
dashed, the zero contour is thicker, and the contour interval is nonuniform. 



Fig. 2. Net radiative heating rates (K d* 1 ) for the nominal IRIS temperatures 
and dust abundances, calculated by summing the solar heating and thermal 
cooii ng rates. Negative contours are dashed, the zero contour is thicker, and the 
contour interval is nonuniform. 
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Elg.3. Mass-weighted stream function, >P m in units of 10 7 kg s- 1 . Positive Rg.5. Eliassen-Palm flax divergence ((pacos (Ji)- 1 V • £)inuoits of m-H d 1. 
values represent clockwise flow, negative contours are dashed, the zero contour Negative contours are dashed, and the zero contour is thicker, 

is thicker, and the contour interval is Donuniform. 
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fig. 4. (a) Diabatic meridional velocity, in m s -1 . Positive values represent 
northward winds, negative contours are dashed, the zero contour is thicker, and 
the contour interval is nonuniform. (b) Diabatic vertical velocity, in cm s~ l . 
Positive values represent upward winds, negative contours are dashed, and the 
zero contour is thicker. 


m 2 s* 1 , two orders of magnitude smaller than the commonly used 
value. The EP flux divergence is found to have thelargest values at 
high latitudes in the winter hemisphere at altitudes above -40 km, 
which coincides with the region of anomalously warm temperatures 
described in [1], 

The diabatic meridional circulation during non-dust-storm con- 
ditions can transport a 1-pm dust particle over global scales, but is 
not intense enough to support 5-fjm particles over a large altitude 
range. The strong meridional circulation also has important impli- 
cations for the role of atmospheric transport in the seasonal cycle of 
water on Mars. At lower atmospheric levels, the meridional circu- 
lation carries water vapor into the equatorial regions, where the 
prevailing winds blow upward. The rising motion leads to ice-cloud 
formation and eventual precipitation of the largest particles, in- 
creasing the water abundance at low latitudes. The water that does 
not precipitate out is transported poleward (as either small ice par- 
ticles or vapor) by the high-altitude meridional winds. This process 
could provide a high-altitude source of water vapor for the formation 
and maintenance of the north polar hood. 

References: [1] Santee M. and Crisp D. (1993) JGR, 98 f 
3261-3279. [2] Crisp D. (1990) JGR, 95 t 14577-14588. 


MARTIAN MESOSCALE CIRCULATIONS INDUCED BY 
VARIATIONS IN SURFACE OPTICAL AND THERMAL 
CHARACTERISTICS : A NUMERICAL STUDY. T.SiilFand 
H. Savijarvi 2 , Finnish Meteorological Institute (FMI), Department 
of Geophysics, P.O. Box 503, 00101 Helsinki, Finland, University 
of Helsinki, Department of Meteorology (DMUH), P.O. Box 4, 
00101 Helsinki, Finland. 

Introduction: According to Mariner 9 and Viking observa- 
tions the surface albedo and the thermal inertia of Mars* surface 
vary substantially, and fairly large gradients also occur [1], The 
Syrtis Major region is bounded in the west by high-albedo cratered 
terrain and in the east by Isidis Planitia, a high-albedo plateau. Sinus 
Meridiani, centered at 0°W 5°S, is almost surrounded by higher 
albedo regions with sharp boundaries, and Acidalia Planitia, be- 
tween 10° and 50° W and north of 35°N, has sharp albedo boundaries 
to the east and west 
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Observational and modeling studies, e.g. f on Australian dry salt 
lake coasts, have shown that discontinuities and gradients in surface 
properties can induce mesoscale circulations [2], We have used a 
version of the DMUH mesoscale model to simulate atmospheric 
circulations induced by variations in the reflectivity and in the 
thermal inertia of the martian surface. 

The Model: A two-dimensional hydrostatic G-coordinatemodel 

is described in detail in [3] and elsewhere. The terrestrial version of 
the model has been used for studies of sea-breezes [c.g.,4] and the 
one-dimensional and two-dimensional Mars versions for summer 
FBL [5] and slope wind [6] studies respectively. The optical depth 
of the atmosphere is spatially and temporally constant and has 
typically been given the value i = 0.4 corresponding to martian clear 
conditions. The soil temperature is predicted with the force-restore 
method, forced by the model’s radiative and turbulent surface 
fluxes. The soil and surface variations are incorporated as profiles 
of topography, albedo, thermal inertia, and surface roughness. Pro- 
files are created off line and given as input to the model; both 
idealized profiles as well as those based on Planetary Data System 
data (see Kg. 1) have been used. The Coriolis parameter is constant 
and in our simulations corresponds to45°N. The season is L,= 135° 
(northern midsummer). Vertical and horizontal mesh extents have 
been *6-8 km and 600-1000 km respectively. 

Circulations Due to Albedo Variations: In our ideal-case 
simulations we have used albedo profiles consisting of two plateaus 
connected with a spline curve. The results show that a reasonable 
albedo gradient of Aa/Ax = 0.1/300 km (see Fig. 1) creates a soil 
temperature differential, which is at its largest in the early afternoon 
and is of the order of 5-10 K, depending on the latitude and season. 

This temperature variation induces an afternoon and evening 
circulation cell superficially resembling the terrestrial sea-breeze 
circulation cells. Clear upward motion of several centimeters per 
second is observed above the low-albedo/high-temperature region 
and sinking motion above the high-albedo/low temperature region. 
Daytime and evening horizontal wind from the high-albcdo region 
to the low-albedo region can be several meters per second and a 
weaker return flow higher up is also observed. In the early evening 
the cell is decoupled from the stable surface layer and stays until 
after midnight A nocturnal flow near the surface opposite to the 
daytime flow (but slightly smaller in magnitude) is also observed. 
The near-surface winds arc fairly light (typically less than 2 m s- 1 ), 
increasing with height and reaching a maximum at a few hundred 
meters of altitude. With no large-scale wind the cell appears to 
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Fig. 2. Surface temperatures during the sol and the surface albedo profile. 
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Fig. 4. Locations of the vertical motion areas at 1930 LT. The contour values 
are in cm s -1 . 

remain — unlike a terrestrial sea-breeze cell — stationary above the 
albedo gradient region. If large-scale wind is present, the cell is 
advected by it; the displacement of the cell can be some hundreds of 
kilometers before collapse. 

As observed in terrestrial sea-breeze situations, increasing pre- 
vailing wind attenuates the vertical motion and the cell strength [7]. 
The cell is weakened by increased optical thickness due to reduced 
short-wave heating. The horizontal wind maximum is approximate- 
ly halved when optical thickness is increased from 0.3 to 1.2; 
vertical flow is also attenuated. 

Conclusions and Future Prospects: Our simulations show 

indications of the occurrence of surface-induced mesoscale circula- 
tions in the martian atmosphere. Especially albedo variations ap- 
pear to generate regular and clear circulation cells with smaller 
effects due to thermal inertia variations. The cell characteristics 
depend on the albedo profile, season, large-scale wind and the 
atmospheric optical thickness. Possible observational verification 
methods are, e.g., surface wind measurements and identification of 
the vcrticalmotionregions by cloud formation from orbital imagery. 
As at some regions onMars the albedo gradients occurprobably due 
to adjacent dust-covered and clear bedrock areas, the horizontal 
winds due to the albedo differentials may play aiole in dust raising 
at those areas. Dust thus raised would modify the circulation itself. 
To study the interaction between the surface dust and the albedo- 
induced circulations, a future avenue of study is inclusion of dust 
lifting and dynamics into the model 

References: [1] KiefferH. H. et al. (1977 )JGR, 82, 4249-4291. 
[2] Physick W. L. and Tapper N. J. (1990 ) Monthly Wea. Rev., 118, 
1029-1042. [3] Alestalo M. and Savijarvi H. (1985) Tellus, 37 A, 
156-162. [4] Savijarvi H. and Alestalo M. (1988) Contrib . Atmos. 
Phys ., 61, 98-104. [5] Savijarvi H. (1992) Contrib. Atmos. Phys., 
64, 219-229. [6] Savijarvi H. and Siili T. (1993) JAS, 50, 77-88. 
[7] Pielke R.A. (1984) Mesoscale Meteorological Modeling, Aca- 
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TWO-DIMENSIONAL MODELING OF THERMAL 
INVERSION LAYERS IN THE MIDDLE ATMOSPHERE 
OF MARS. B. Theodore and E. Chassefiere, Serviced* Aeronomic 
du CNRS, BP 3, 91371 Verrieres-le-Buisson Cedex, France. 

There is some evidence that the thermal structure of the martian 
middle atmosphere may be altered in a significant way by the 
general circulation motions. Indeed, while it is well known that the 


circulation in the meridional plane is responsible for the reversal of 
the latitudinal thermal gradient at the solstice through the adiabatic 
heating due to sinking motions above the winter pole, here we want 
to emphasize that a likely by-product effect could be the formation 
of warm layers, mainly located in the winter hemisphere, and 
exhibiting an inversion of the vertical thermal gradient 

Observations: A noteworthy result was obtained from hetero- 
dyne measurements using the IRAM 30-m telescope: A mapping of 
the CO(J = 2-1) line was performed during the 1988 opposition, 
while Mars was under northern winter solstice conditions. It al- 
lowed the retrieval of the vertical thermal profile as a function of the 
position on the martian disk, with a spatial resolution of one half 
hemisphere. The winter profiles definitely exhibit a bump near 50 
km altitude, i.e., a relative warming of about 5 K compared to the 
southern profiles. This may be interpreted as the presence of a warm 
layer above the high winter latitudes, implying an inversion of the 
vertical thermal gradient of about 30 K [lj. Such inversions were 
inferred, but with a weaker amplitude, from the presence at the 
equinox and near the equator of water ice clouds detected by solar 
occultations performed aboard the Phobos 2 spacecraft [2]. 

It must be pointed out that this phenomenon is distinct from the 
polar warming reported by Martin and Kieffer [3] from the Viking 
lander data, which is believed to be quite similar to the terrestrial 
sudden stratospheric warming (SSW) and may be linked to the 
planetary waves activity [4] . The warming we observed is probably 
analogous to the terrestrial mesospheric inversions detected below 
the mesopause by lidar [5] or from space [6], and it was shown [5] 
that their presence is a response to the dynamical forcing provided 
by the breaking of internal gravity waves. 

It is well accepted that gravity waves play a key role in atmos- 
pheric dynamics. The energy budget of the atmosphere requires the 
presence of a strong dynamical forcing in order to be balanced. The 
breaking of gravity waves could play this role: the momentum 
deposition drags the flow and generates turbulence. The former 
closes the zonal jets near the mesopause and leads to the formation 
(or the persistence) of a meridional circulation cell because of the 
Coriolis torque, while the latter leads to a local turbulent heating and 
creates a jump in the eddy diffusion coefficient vertical profile that 
seems to be realistic [7]. Yet the main effect on the thermal structure 
of the atmosphere is due to vertical motions in the meridional cell 
that produce adiabatic warming or cooling; it is how the latitudinal 
temperature gradient is inverted. This phenomenon is well known 
in the terrestrial mesosphere and it seems likely in the case of Mars. 
Indeed, Deming et al. [8] observed that the middle atmosphere is far 
from radiative equilibrium above polar regions, unlike equatorial 
ones. The formation of inversion of the vertical temperature gradi- 
ent is believed to follow the same mechanism, but implies the 
momentum deposition to occur in a narrow layer instead of the 
whole atmosphere, leading to local heating. 

Modeling: In order to assess the relative importance of the 

various processes involved in this phenomenon, we adapted a two- 
dimensional model of the terrestrial mesosphere, initially devoted 
to the study of gravity waves-mean flow interactions in the middle 
atmosphere of the Earth and developed by Hauchecome and Maillard 
[9]. We improved this code by the addition of the vertical and 
meridional advection terms and by optimizing the filtering [10]. 
This work may be considered as following the exploratory work of 
Barnes [1 1], who showed that gravity wave activity could lead to a 
strong warming in the middle martian atmosphere. 
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Our code solves the zonally averaged primitive equations in 
spherical coordinates. Gravity wave effects appear in terms of 
friction in the zonal momentum equation (the waves are assumed to 
be purely zonally propagating), turbulent heating in the thermody- 
namic equation, and vertical diffusion of all the components. Mo- 
lecular diffusion is neglected and filtering plays the role of the 
horizontal diffusion. Radiative processes are parametrized with a 
Newtonian cooling that relaxes permanently the thermal field to- 
ward the initial “radiative” field. 

The gravity wave effect parameterization is basically the same as 
formulated by Lindzen [12], using the saturation hypothesis: waves 
propagate upward until they reach a level where they become 
unstable. Above this level, up to the critical level (determined as the 
level where the wave phase speed is equal to mean flow velocity), 
the wave amplitude is assumed to remain constant; just enough 
turbulence is produced to prevent the growth of the wave. We 
assume the wave to be not damped by the radiative processes or by 
the eventual background turbulence below the breaking level. Holton 
and Zhu [13] and Barnes [11] showed that this has little influence, 


but it is perhaps a crude assumption. In order to simulate a spectrum 
of waves we use three phase velocities associated with three wave 
amplitudes, and we make the assumption that the wave generation 
is latitudinally uniform. 

As the program runs, the breaking of gravity waves leads to the 
closure of the zonal wind jets through the drag exerted near the 
mesopause (Fig. 1). The strong meridional winds induce a deep 
penetration of the summer jet in the northern hemisphere (up to 
20°N), as suggested by wind measurements in the middle martian 
atmosphere [14]. Furthermore, a warm layer develops above the 
winter polar regions, leading to an inversion of the vertical thermal 
gradient of about 30 K. This inversion occurs by adiabatic heating 
in the sinking branch of the meridional cell. A noteworthy result was 
obtained: forsome wave characteristics, a wanning occurs above the 
equatorial regions. Its amplitude depends on several parameters, 
like the wave characteristics as well as the value of the Newtonian 
cooling coefficient (which may be related to the optical depth of the 
atmosphere). It is due to the action of gravity waves on the tongue 
of retrograde winds entering the winter hemisphere; indeed, they 





Fig. I. A latitude-altitude cross section of the zonal, meridional, and vertical wind and temperature. 
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tend to counteract the meridional circulation locally. The amplitude 
of the inversions strongly depends on the wave parametere and they 
can vanish if the wave amplitude is too weak. 

Conclusions: This work reinforces the belief in favor of a 
strong wave activity in the middle atmosphere of Mars at planetary 
scales. Using a simple numerical simulation of the interaction 
between internal gravity waves and general circulation, which must 
be considered at the present stage as a heuristic model dedicated to 
the understanding of new phenomena, still poorly observed and 
constrained, it has been possible to successfully interpret global- 
scale millimeter-wave observations of wind velocities and thermal 
structure in the martian middle atmosphere. The strong retrograde 
winds observed in the whole s umm er hemisphere up to the low 
latitudes of the winter hemisphere [14], and the relative atmospher- 
ic warming above the winter hemisphere near 50 km can be simul- 
taneously explained by the effect of breaking gravity waves. 

It is noteworthy that global- scale observations of temperature 
and wind fields in the middle atmosphere of Mars may be explained 
by using our schematic two-dimensional model. In spite of the crude 
initialization and other simplifying assumptions relative to the 
modeling of wave propagation and breaking, a consistent picture 
may be found. This shows that, as originally suspected by Barnes 
[11], gravity waves could play a major role in the setting of the 
dynamical regime of the martian middle atmosphere. In addition, 
new phenomena, like equatorial temperature inversions, are pre- 
dicted. It must be noted that such inversions are observed in the 
Earth*s mesosphere [6], but still not understood, suggesting that the 
study of this phenomenon might help modelers of the Earth’s meso- 
sphere. The martian atmosphere, simpler than the terrestrial atmo- 
sphere and expected to be dynamically and thermally more con- 
trasted, could turn out to be a good natural laboratory for the study 
of wave activity. 

New millimeter-wave observations, obtained in December 1992, 
will provide additional information on the dynamics of the middle 
atmosphere of Mars near the equinox. The Mars Observer mission 
will bring substantial information on the thermal structure of the 
middle martian atmosphere at all seasons, mainly through PMIRR 
and Radio Science experiments. The use of the present model for the 
interpretation of PMIRR data will constitute a unique opportunity 
to test in a detailed way and to improve it, for example, by introduc- 
ing more realistic initialization of the wind and thermal fields. 
Implications for the vertical distribution of some photochemically 
active constituents like ozone will be tested by using the results of 
experiments mounted onboard Mars 94/96 missions, like the 
SPICAM experiment 
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THREE-DIMENSIONAL NUMERICAL SIMULATION OF 
THERMAL TIDES IN THE MARTIAN ATMOSPHERE. 
R. J. Wilson, Geophysical Fluid Dynamics Laboratory/NOAA, 
Princeton University, P.O. Box 308, Princeton NJ 08542, USA 

Analysis of the Mariner 9 IRIS and Viking Orbiter IRTM tem- 
perature data and of the time series of surface pressure from the 
Viking landers has indicated that atmospheric thermal tides are a 
significant component of the martian general circulation. Classical 
tidal theory, considering only the Sun synchronous components, has 
been shown to at least roughly account for the amplitudes of the 
observed diurnal and semidiurnal surface pressure oscillations. In 
particular, Zurek and Leovy [1] have demonstrated that classical 
tidal theory can reproduce the observed dependence of the strength 
of the semidiurnal tide on aerosol optical depth. The high-amplitude 
topography of Mars, however, can cause longitudinal distortions in 
the thermo tidal forcing, leading to additional tidal modes that are 
not Sun synchronous. In particular, observation [2] and theory [3] 
have suggested the likely presence of a diurnal, zonal wavenumber- 
1 K wave. Hamilton and Garcia [4] noted that the wavenumber-1 K 
normal mode has a period of close to 24 hr and may be detectable 
in surface pressure observations. Zurek [5] has discussed the possi- 
bility of this mode being resonantly excited for various atmospheric 
temperature structures. The semidiurnal wavenumber-2 K mode, 
with roughly a 12-hr period, could also be resonantly excited. A 
numerical model with finite amplitude topography and a self-con- 
sistent thermotidal forcing resulting from the daily heat flux at the 
surface and atmospheric absorption of solar radiation by aerosols is 
used here to explore the influence of topography on martian thermal 
tides and to examine the possibility of Kelvin wave resonance. 

A fully nonlinear, three-dimensional, global primitive equation 
model has been constructed to address some of the limitations of 
linear theory. In order to obtain realistic thermotidal forcing, the 
model physics include a three-level soil model, a Monin-Obukhov 
surface drag formulation, Richardson number-dependent vertical 
diffusion, and convective adjustment of heat (but not momentum). 
The model does not yet provide for the C0 2 condensation circula- 
tion. Experiments with the NASA Ames martian GCM [6] suggest 
that this omission is not important in determining the zonal mean 
wind and temperature fields except at the surface in the vicinity of 
the polar caps and, hence, is not of great significance for tidal 
studies. Atmospheric heating by solar radiation is due to absorption 
in the NIR bands of C0 2 and by aerosols. Heating by aerosols is 
calculated using a delta-Eddington approximation for aerosol scat- 
tering. The optical parameters are those suggested by Pollack et al. 
[7]. Cooling by dust at longwave frequencies is accounted for in the 
same fashion as in Haberle et al. [8] using an emissivity formulation 
for the region outside the 15-pm band. The IR effects of C0 2 are 
based on the simple 15-pm C0 2 band model of Goody and Belton 
[9]. At present the model is usually run with a uniform surface 
thermal inertia and albedo except where the surface temperatures 
imply the formation of C0 2 ice. The model has been run with 25 or 
40 levels (typically 40) in the vertical, with the domain extending 
to roughly 45 or 90 km in height To date, all runs have used a 5° x 
6° horizontal resolution. 

In each of the experiments described here, the aerosol has a 
horizontally uniform distribution with visible optical depth chosen 
to be 0.3, 1.0, or 2.0. The zonally averaged winds and temperatures 
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agree quite well with those from the NASA Ames GCM [6]. In 
particular, the circulations have a similar dependence on the season 
and on dust optical depth. The two models exhibit comparable 
amplitudes of diurnal temperature variability up to 40 km. 

The first experiments employed an idealized wavenumber-2 
topography with a cosine dependence in latitude. A series of simu- 
lations was carried out, varying the topographic amplitude from 
03-3 km and using different dust optical depths. Figure 1 displays 
the results of a harmonic analysis of the model surface pressure for 
2-km amplitude topography. The eastward and westward propagat- 
ing components of the zonal wavenumber-1 diurnal tide are shown, 
as well as the corresponding wavenumber-2 components of the 
semidiurnal tide. The figure clearly shows the presence of an 
eastward propagating wavenumber- 1 K wave with diurnal fre- 
quency. The strength of this Kelvin wave response varies approxi- 
mately linearly with topographic amplitude over the range of ampli- 
tudes considered, and is fairly insensitive to atmospheric tempera- 


ture variations associated with variable dust loading. As expected, 
the westward (Sun-synchronous) tides show a pronounced sensitiv- 
ity to increased dust, with the semidiurnal tide having relatively 
greater sensitivity. The strength of the semidiurnal wavenumber- 
2 K wave is apparently highly dependent on atmospheric thermal 
structure as well as topographic amplitude. With only wavenumber - 
2 topography present, this presumably resonant mode must have 
been forced by nonlinear interactions. 

A latitude-height plot of the amplitude of the wavenumber-1, 
eastward-propagating, diurnal component of zonal velocity has a 
classical Kelvin wave signature in the lower atmosphere, but dis- 
plays a pronounced tilt with height towards the westerly jet in the 
winter hemisphere. 

Simulations with highly smoothed versions of realistic topogra- 
phy bear out the results of the idealized topography cases. Again, the 
zonal wavenumber-2 component is crucial for the Kelvin wave 
excitation. Figure 2 shows the spatial variability of the amplitude of 
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Fig. 1. Harmonic analysis of the surface oscillations for the case of 2.0- km idealized topography at the time of southern hemisphere solstice. 
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Fig. 2. Diurnal and semidiurnal pressure amplitudes as simulated with 
smoothed realistic topography at the time of southern hemisphere solstice. Dust 
optical depth is 0.3; amplitude is in microbars. 


the diurnal and semidiurnal surface-pressure oscillations for a case 
with dust optical depth of 0.3. There is a considerable degree of 
topographic modulation in the diurnal tide and a fair degree in the 


VERTICAL TRANSPORT OF WATER IN THE MARTIAN 
BOUNDARY LAYER. A. P. Zent*. R. M. Haberle 2 , and H. C. 
HoubenVSETI Institute andNASA Ames Research Center, Moffett 
Field CA 94035, USA, 2 NASA Ames Research Center, Moffett 
FieldCA94035, USA, 3 SpacePhysics Researchlnstitute andNASA 
Ames Research Center, Moffett Field CA 94035, USA. 


We are continuing our examination of the transport of H 2 0 
through the martian boundary layer, and we have written a one- 
dimensional numerical model of the exchange of H 2 0 between the 
atmosphere and subsurface of Mars through the planetary boundary 
layer (PBL) [1], Our goal is to explore the mechanisms of H 2 0 
exchange, and to elucidate the role played by theregolith in the local 
H 2 0 budget The atmospheric model includes effects of Coriolis, 
pressure gradient and frictional forces for momentum, as well as 
radiation, sensible heat flux, and advection for heat The model 
differs from Flasar and Goody by use of appropriate Viking-based 
physical constants and inclusion of the radiative effects of atmo- 
spheric dust We specify the pressure gradient force or compute it 
from a simple slope model. The subsurface model accounts for 
conduction of beat and diffusion of H 2 0 through a porous adsorbing 
medium in response to diurnal forcing. The model is initialized with 
depth-independent H 2 0 concentrations (2 kg nr 3 ) in the regolith, 
and a dry atmosphere. The model terminates when the atmospheric 
H 2 0 column abundance stabilizes to 0.1% per soL 

Results suggest that in most cases the flux through the martian 
surface reverses twice in the course of each sol. In the midmorning, 
the regolith begins to release H 2 0 to the atmosphere, and continues 
to do so until midafternoon, when it once more becomes a sink. It 
remains an H 2 0 sink throughout the martian night In the early 
morning and late afternoon, while the atmosphere is convective, the 
atmosphere supplies H 2 0 to the ground at a rapid rate, occasionally 
resulting in strong pulses of H 2 0 into the ground. The model also 
predicts that for typical conditions, perhaps 15-20 sols are required 
for the regolith to supply an initially dry atmosphere with its 
equilibrium load, although more than half of the equilibrium load 
typically appears in the atmosphere on the first soL 

The effects of surface albedo, thermal inertia, solar declination. 


semidiurnal tide, which clearly is a departure from the longitudi- 
nally uniform amplitudes of classical tidal theory. An experiment 
with wavenumber-2 variation in the thermotidal forcing due to 
surface albedo variation (0.15-0.35) produced only a rather weak 
Kelvin wave response relative to that with topography. A series of 
simulations with a varying gas constant (which should change the 
resonant frequency of the Kelvin mode) su ggests that the solar tides 
in the martian atmosphere are indeed within a resonance regime. 

Work is in progress to examine the structure of linear tidal modes 
in the presence of strongly varying basic states and to determine the 
effects of the diurnal cycle on the simulated general circulation. 

References: [1] Zurek R. W. and Leovy C. B. (198 1) Science, 

213, 437-439. [2] ConrathB. J. (1976)/AS, 33, 2430-2439. [3] Zurek 
R. W. (1976) JAS, 33, 321-337. [4] Hamilton K. and Garcia R. R. 
(1986) JGR, 91, 11867-11875. [5] Zurek R. W. (1988) JGR, 93, 
9452-9462. [6] Habcrle R.M. et al. ( 1993) JGR, 98, 3093-3123. [7] 
PoUack J. B. et al. (1979) JGR, 84, 2929-2946. [8] Haberle R. M. 
ct al. (1982) Icarus, 50, 322-367. [9] Goody R. M. andBelton M. 
J. S. (1967) Planet Space ScL, 15, 247-256. 


atmospheric optical depth, and regolith pore structure have been 
explored. Increased albedo cools the regolith, so less H 2 0 appears 
in the atmospheric column above a bright surface. The friction 
velocity is higher above a dark surface, so there is more diurnal H 2 0 
exchange; relative humidities are much higher above a bright sur- 
face. Thermal inertia (I) affects the propagation of energy through 
a periodically heated homogeneous surface. Our results suggest that 
higher thermal inertia forces more H 2 0 into the atmosphere because 
the regolith is wanner at depth. Surface stresses are higher above a 
low I surface, but there is less diurnal exchange because the atmos- 
phere is dry. The latitude experiment predicts that the total diurnal 
insolation is more important to the adsorptively controlled H 2 0 
column abundance than the peak daytime surface temperature. Fogs 
and high relative humidity will be far more prevalent in the winter 
hemisphere. The dust opacity of the atmosphere plays a very signifi- 
cant role; the PBL height, column abundances, relative humidity, 
and surface stresses all increase very strongly as the optical depth 
approaches zero. The dust opacity of the atmosphere must be 
considered in subsequent PBL models. 
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Current plans for upgrades include making the thermal proper- 
ties of the regolith temperature dependent, and including the latent 
heat of phase change for atmospheric and subsurface H 2 0. While all 
previous simulations terminated after equilibration between the 
atmosphere and subsurface, we are currently allowing the model to 
time-march through several degrees of L* in order to track seasonal 
changes in the vertical distribution of atmospheric and subsurface 
H 2 0. We are also exploring the effects of changes in atmospheric 
opacity. The equilibrium column abundance of H 2 0 in a clear 
atmosphere is much higher than in a dusty atmosphere. However, it 
is notclear that an increase in optical depth, such as due to the arrival 
of a dust storm, could remove water from the atmosphere, since the 
diurnal boundary layer collapses as the optical depth increases. This 
may be analogous to the late afternoon collapse of the boundary layer 
due to increased static stability, which effectively strands a consid- 
erable amountof water at altitude, even though it is out of adsorptive 
equilibrium with the cold surface. 

References: [1] Zent A: P. et aL (1993) JGR f 98, 3319-3338. 


THE MARS AERIAL PLATFORM (MAP) MISSION. 
R. Zubrin 1 , S. Price 1 , T. Gamber 1 , B. Clark 1 , R. Haberle 2 , and J. 
Cantrell 3 , ! Martin Marietta Astronautics, P.O. Box 179, Denver 
CO, 80201, USA, *NASA Ames Research Center, Moffett Held 
CA, 94035, USA, 3 Space Dynamics Lab, Utah State University, 
Logan UT 84322, USA. 

The Mars Aerial Platform (MAP) mission is a conceptual design 
for a low-cost, Discovery-class mission whose purpose would be to 
generate tens of thousands of veiy-high-resolution (20 cm/pixel) 
pictures of the martian surface, map the global circulation of the 
martian atmosphere, and examine the surface and subsurface with 
ground-penetrating radar, infrared spectroscopy, neutron spectros- 
copy, and other remote sensing techniques. The data would be 
acquired by instruments that are carried by balloons flying at a 
nominal altitude of about 7 km over the martian surface. Because 
new balloon and microspacecraft technology is now available, the 
balloon probes could be quite long-lived, lasting hundreds or even 
thousands of days, producing an immense science harvest in the 
process. Together with the Mars Environmental Survey (MESUR) 
surface network science mission, MAP would revolutionize our 
knowledge of the Red Planet 

Mission Description: The MAP mission will be carried out as 
follows: A low-cost launch vehicle such as a Delta 7925 (1000 kg 
TMI) is used to propel a small spacecraft carrying eight entry 
capsules onto trans-Mars injection. Approaching Mars, a spinner is 
used 10 days prior to arrival to release the capsules outward so that 
they enter Mars’ atmosphere at widely dispersed locations. Each 
capsule then enters the atmosphere, and deploys a parachute that 
slows it down to the point where a balloon can be inflated. The 
inflation is accomplished during descent, so that no landing system 
is required (the practicality of accomplishing this has been demon- 
strated on Earth at altitudes of up to 150,000 ft, where atmospheric 
density is similar to Mars). After the balloon is inflated, the para- 
chute, capsule, and inflation equipment is jettisoned, and used to 
land a meteorology package consisting of pressure and temperature 
sensors, battery, and transmitter on the surface. Each of the balloon 
probes will then commence their float around Mars at an altitude of 
about 7 km (23,000 ft) above the mean surface level. The probes, 


which will never land in the course of their long-duration aerial 
cruise, will cany a gondola payload of 8 kg, which includes 2 
cameras, 1 kg of atmospheric science and additional instruments, 
data recording and transmitting equipment, a rechargeable battery, 
and solar array. 

When it is daylight two pictures are taken simultaneously every 
15 minutes. One is taken with a high-resolution CCD black and 
white camera with a nominal resolution of 20 cm per pixel; with a 
1024 x 1024 pixel field, this gives a field of view 204 m on a side. 
The other is taken with a color camera at a moderate resolution of 
10 m per pixel; this gives a field of view 10.2 km on a side. The two 
cameras are aligned so that the high-resolution image is located at 
the center of the medium-resolution picture, whose features in turn 
can be used to locate the region studied on a map of the planet The 
cameras are mounted on a vertical pan, so that either vertical, 
horizontal, or oblique pictures can be taken. The data is compressed 
and stored and then is periodically uplinked to an orbiting satellite, 
which could either be Mars Observer or the MESUR comsat, and is 
then transmitted back to Earth. Under typical conditions, each 
aerostat will be able to transmit 40 picture pairs per day. It is 
anticipated that the science return of this mission will be large, with 
32,000 high-resolution pictures returned to Earth for each 100 days 
of operation of an eight balloon fleet, plus an equal number of 
pictures at resolutions superior to the best Viking images. The high- 
resolution pictures will have a factor of 7 greaterresolution than the 
best pictures returned from Mars Observer. Additional science 
return will result from use of other instruments, such as a ground- 
penetrating electromagnetic sounder. Tracking the balloons will 
yield knowledge of Mars* global atmospheric circulation. The cost 
of the mission is anticipated to be low, falling below the $150 
million (without launch) ceiling of the Discovery-class guidelines. 

Science Return: Some of the science areas that will be ben- 

efited by the MAP mission will include the following: 

1. Geologic science. The greatest benefit from high-resolu- 
tion images is the information that could be gained on surface 
processes. In contrast to Mars Observer, whose highest-resolution 
images will show objects as small as a mid-sized car, the cameras 
on the MAP gondola will be able to reveal surface features the size 
of a cat. Such images will give us submeter surface morphology, 
block size distributions, and surface textures, all of which provide 
insight into volcanic processes, sedimentary deposits, fluvial histo- 
ry, and aspects of impact cratering. In addition, colian features such 
as ripple fields, individual dunes, small channels, etc., will be 
visible, providing direct evidence for the link between the atmo- 
sphere and the surface. The scale of some of these is below Mars 
Observer resolution. Assessment of these and other features that are 
submeter in size will provide a whole new suite of information on 
surface/atmosphere interactions. 

2. Atmospheric science. The MAP mission provides a unique 
in situ platform for atmospheric science. Each balloon will float at 
a fixed density level (0.008 kg nr 3 * or about 7 km above the mean 
surface) for at least 100 days. If the balloons were to average 
10 m s- 1 , aconservative value, then thousands of kilometers will be 
traversed. Depending on the circulation, it is possible that some of 
the balloons will circumnavigate the planet several times. Thus 
global sampling of the atmosphere is possible. 

Measurements that can be taken aboard the MAP balloons 
include temperature, pressure, horizontal wind speed and direction 
(via tracking), dust size and composition, and water vapor concen- 
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tration. Atmospheric science goals that could then be addressed 
include the critically important surface boundary' layer and its 
diurnal variability, atmospheric dynamics, circulation systematic^, 
aerosol composition and properties, and the seasonal cycles of dust 
and water. Tracking of the balloons would give information on 
winds and circulation systems that could be used to test the compu- 
tational global circulation models currently being developed by 
NASA. Slope wind systems may be particularly important on Mars 
and the MAP balloons should be able to help characterize their 
structure and dynamics. Aerosols, dust in particular, are known to 
be ubiquitous in the martian atmosphere, yet their composition, 
size, and shape are uncertain. MAP balloons could measure these 
properties directly and thereby greatly increase our understanding 
of these climatically significant particles. Water vapor concentra- 
tions could be measured, furthering our understanding of the nature 
and distribution of near surface water reservoirs. 

3. Exploratory science . Infrared detectors on the balloons 
would be able to detect local geothermally heated hot-spots on the 
martian surface, should any exist and fall within several kilometers 
of the probe’s ground track. Neutronic devices would be able to 
detect local concentrations of ice, hydrated minerals, or deposits of 
carbonates. Ice layers 10 m or more below the surface may also be 
revealed by electromagnetic sounding using the radar altimeter, 
which can function as a simple ground-penetrating radar. 

The discovery of such features would be of very high value for 
exobiological investigations. 

4. Development of Mars engineering models. Numerous U.S . 
and international groups have a high interest in developing engi- 
neering models of the martian surface. A recurring theme is the need 
for imaging of the surface at the scale that could be provided by 
MAP. In particular, knowledge of size-block distributions, slopes, 
surface discontinuities, etc., at a scale important to landers, roving 
vehicles, and other landed mission operations is of high interest 
The advantage of MAP is that a wide variety of terrains could be 
sampled to provide input to engineering models. The Mars Observer 
camera has a resolution of about 1.5 m/pixel, allowing it to resolve 
features about 4 m across. Such resolution is too gross for use in 
assessment or avoidance of landing hazards facing Mars surface 
science missions . S ubmeter scale imaging is also needed to provide 
context for the activities of a landed rover mission. This context 
would permit significantly richer surface operations activities, and 
a concurrent improvement in the caliber of science conducted. The 
results from the MAP aerostats will be obtained concmrently with 
data returned by the meteorology packages delivered to the surface 
by the entry probe parachutes. This will enhance the potential 
science return by the synergy that will result from concurrent 
measurements from the surface network with the balloon drift 
measurements. This will provide more information on the atmo- 
spheric structure than could be obtained by either a balloon or 
surface network acting alone. The large array of potential science 
instruments that could be carried by MAP are s umm arized in 
Table 1 . There are far too many to be all carried in the same gondola; 
however, if the probes sent to Mars arc divided into four “science 
sets” (which we label “A,B,C, D”) with two probes per set for a 
Delta launch, then the full range of potential investigations could be 
carried out. To minimize power consumption, all instruments would 
be cycled on and off in coordination with the camera. 

Balloon Flight Profile: At the time of entry, the balloons may 

be targeted to widely dispersed locations, typically spread over 130° 


of latitude. The balloons will then commence to float at the velocity 
of the martian winds at their altitude. Current models suggest that 
the wind will cany the balloons primarily in the west-east direction 
with typical velocities on the order of 30— 60 km hi -1 . At such speeds 
each aerostat could be expected to circumnavigate the planet every 
10-30 days. There will also be some north-south migration, which 
will probably cause some balloons to wander into and out of Man' 
polar regions during portions of the mission. This may afford an 
opportunity to examine regions of the planet that cannot be targeted 
for direct probe entry during a given launch opportunity. If an 
aerostat should wander into a polar region during the Arctic night, 
gondola operations will be curtailed until it leaves. Characteristic 
aerial velocities of the balloons suggest that such a shutdown should 
last no more than a few days at most. Since the thermal protection 
power needed by the MAP gondola is provided by three Radioiso- 
tope Heating Units (RHUs), which require no electric power, the 
gondola is survivable even if solar incidence should drop far below 
anticipated levels for an extended period. 

The balloons have a nominal float altitude during the northern 
summer of 7 km, with a diurnal altitude variation of a few hundred 
meters. However, because of the evaporation of C0 2 from Mars* 
south polar cap, the average atmosphere during the course of the 
year is thicker than that of the northern summer, and so the balloons 
on the average will actually float at atypical altitude of about8.5 km, 
clearing over 97% of all martian terrain. Having such a large fraction 
of the planet safe for travel gives each balloon fair odds for a long- 
duration flight Although an eventual crash is certainly possible, 
chances for balloon loss through collision with obstacles is dimin- 
ished by the tendency of balloons to follow the streamlines of air 
which must, of necessity, go around all obstacles the airflow encoun- 
ters. 

Probe Navigation: Latitude of the MAP aerostats can be 
determined by the time of comsat overpass, which will occur 2-5 
times each sol. In addition, the perpendicular distance of the probe 
from the comsat groundtrack can be determined by the Doppler 


TABLE 1. Candidate MAP instruments. 


Instrument 

Mass 

(kg) 

Power 

(W-h/sol) 

Data 

(Mbar/sol) 

Science Set 
(example) 

Camera (BW or 3-color) 

0.6 

10 

26.0 

AB,C,D 

Temperature and pressure 

0.1 

0.01 

0.02 

AB, C, D 

Sun sensor 

0.1 

0.01 

0.02 

A B, C,D 

Accelerometer 

0.1 

0.1 

0.2 

A C 

Neutron spectrometer 

0.2 

0.03 

0.06 

A, B,C 

Gamma ray spectrometer 

1.0 

2.0 

5.0 

B 

Electromagnetic sounder 

0.4 

8.0 

1.0 

A C 

Microwave radiometer 

0.3 

0.1 

0.02 

A 

Magnetometer 

0.2 

2.0 

0.02 

A D 

Infrared radiometer 

0.1 

0.01 

0.02 

A C,D 

Infrared spectrometer 

0.5 

2.0 

5.0 

C 

Nephelometer 

0.2 

0.5 

1.0 

D 

Hygrometer 

0.2 

0.02 

0.02 

D 

Dust composition analyzer 

0.4 

0.1 

0.02 

D 
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shifting of the signal the comsat receives from the Probe. In addition, 
longitude can be determined twice a day by noting the time of 
sunrise and sunset If longitude is known, then a Sun sensor can be 
used to determine latitude and gondola orientation. The aerostat 
location can also be determined by comparing its 10 x 10 -km 
medium-resolution images with Mars Observer and Viking surface 
imagery. Since the position of the probe is periodically known 
precisely, the position of the probe at all times can be known if 
overlap can be maintained between the areas captured by the 
medium-resolution images. Since the probe can know its own 
location using such techniques, probes can be preprogrammed to 
take horizontal photographs when they find themselves in a certain 
region and pointing (the balloons are expected to slowly rotate) in 
a given direction. Thus it will be possible to obtain horizontal or 
oblique panoramic images of many of the major sights on Mars, such 
as Olympus Mons or the Valles Marineris. 

Conclusion: We conclude that the MAP mission offers the 

potential for a very high science return mission at modest cost, and 
recommend that appropriate steps be taken to implement it at an 
early date. 

MARS OBSERVER: APPLICATIONS TO ATMOS- 
PHERIC TRANSPORT. R. W. Zurek and D. J. McCleese, Jet 
Propulsion Laboratory, Califomialnstitute of Technology, Pasadena 
CA 91109, USA. 

The Pressure Modulator Infrared Radiometer (PMERR) is an 
atmospheric sounder designed to observe temporal and spatial 
variations of water vapor and of dust suspended in the Mars atmos- 
phere, to characterize the planetary-scale thermal structure and 
circulation of the atmosphere, and to quantify the polar radiative 
balance [1]. These measurements are fundamental to understanding 
the seasonal cycles of dust, of water, and of C0 2 on Mars and, in 
particular, to determining the role of atmospheric transport Using 


measurements in eight nairow-band infrared spectral regions and 
one broadband visible channel, the PMIRR investigation teams at 
JPL and Oxford University will derive vertical profiles of atmos- 
pheric temperature, extinction due to suspended dust and water 
vapor concentration, as well as locations of CO z and H 2 0 ice clouds. 
These data will be used in a variety of ways to address issues of 
atmospheric dynamics and transport 

Three topics will be emphasized here: ( 1) the expected precision 
of the retrieved profiles of temperature, dust extinction, and water 
vapor, including plans for validating the profiles; (2) the observa- 
tional strategy, designed to best use PMIRR' s two-axis scan mirror, 
as deployed in thcMars Observer mapping orbit; and (3) approaches 
to mapping the atmospheric fields globally and the derivation of key 
meteorological fields related to estimating atmospheric transport 

The underlying philosophy of the observing strategy is to map the 
atmospheric fields as continuously and as systematically as pos- 
sible, since much of what we want to learn comes from observing 
variations with both space and time. Changes in the mapping 
strategy as a function of season (with a planet-encircling dust stoim 
counted as a “seasonal” event) will be discussed. The individual 
vertical profiles will be mapped using different algori thms to pro- 
vide synoptic global fields of temperature, of dust, and, where 
possible, of water vapor . These mapped fields can be used to 
directly derive estimates of wind shear and of radiative forcing; they 
and/or the retrieved vertical profiles will also be provided to several 
groups for assimilation of the retrieved fields into data-model 
blends that will provide estimates of the meteorological fields, 
including fields that cannot be derived directly (e.g., reference-level 
winds and vertical velocities). While such data ass imil ation tech- 
niques can, in principle, provide “optimal” estimates of the fields, 
there are many issues yet to be resolved, including proper character- 
ization of observational errors and of deficiencies in the dynamical 
models. 

References: [1] McCleese D. J. et al. (1992) JGR , 97, 

7735-7757. 
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